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Abstract. The notion of r-crossing and r-nesting of a complete match- 
ing was introduced and a symmetry property was proved by Chen, Deng, 
Du, Stanley and Yan in 2007. We consider random matchings of large 
size and study the maximal crossing and the maximal nesting. It is 
known that the marginal distribution of each of them converges to the 
GOE Tracy- Widom distribution. We show that the maximal crossing 
and the maximal nesting becomes independent asymptotically, and eval- 
uate the joint distribution for the Poissonized random matchings explic- 
itly to the first correction term. This leads to an evaluation of the 
asymptotic of the covariance. Furthermore, we compute the explicit 
second correction term in the distribution function of two objects: (a) 
the length of the longest increasing subsequence of Poissonized random 
permutation and (b) the maximal crossing, and hence also the maximal 
nesting, of Poissonized random matching. 



1. Introduction 

Let A4n be the set of complete matchings of [2n]. The size of is (2n — 
1)!!. It is well-known that the number of complete matchings of [2n] with no 
crossings equals the nth Catalan number C„, as is the number of complete 
matchings with no nestings. In |15j . a notation of r-crossing and r-nesting 
was introduced: given a complete matching M = {(ii, ji), . . . , (i„, j„)} G 
■Mn, {iisi,jsi),---,iisr,jsr)} IS callcd an r-crossing if i^^^ < ig^ < ■ ■ ■ < 
isr < jsi < ■ ■ < jsr and an r-nesting if i^i < isa < ' ' ' < ^sr < Jsr < ' ' ' < 
js2 < jsi- Let cr„(M) be the largest number k such that M has a /c-crossing 
(maximal crossing) and nen{M) denote the largest number j such that M 
has a j-nesting (maximal nesting). See Figure [T] for an example. Various 
combinatorial properties of cr„ and ne^ were studied by Chen, Deng, Stanley, 
and Yan in |15j . This paper subsequently generated a flurry of research 
concerning crossings and nestings of many combinatorial objects: see, for 
example, |42] and also the survey article |48| . 

We may equip Mn with the uniform probability and regard cr„, and ne„ 
as random variables. Let M he a Poisson random variable with parameter 
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1 2 3 4 5 6 7 S y 10 11 12 



Figure 1. A complete matching M of [12]. In this sam- 
ple CTeiM) = 4, achieved by {(1, 6), (2, 7), (4, 9), (5, 10)}, and 
nee{M) = 2, achieved by {(3, 11), (4, 9)}. 



and consider matchings of random size distributed as 2Af. Let CRt and 
NEj denote cr_^ and nej\f, respectively. The object of this paper is to study 
the asymptotics of CRj and NEt as t — )• oo. 

One of the main results of [12] is that the joint distribution of cr„ and ne,i 
are symmetric. Hence CRj and NE^ are symmetrically distributed. The limit 
of the marginal distribution of NE^ can be obtained by noting a bijection be- 
tween matchings and fixed-point-free involutions. Let Inv„ be the set of per- 
mutations of size 2n consisting of only 2-cycles. To cj € Inv„ whose cycles are 
(ii, ji), . . . , (injn), associate the complete matching {(hji), (injn)}- 
This gives a natural bijection ip from Invn onto A4n- Moreover, if we define 
in{cr) as the length of the longest decreasing subsequence of o" G Inv„, it 
is easy to check that in{cr)/2 = ne„((/j((T)). The limiting distribution of £„, 
and also of ij\f were obtained obtained earlier in [U |9] . From this and the 
symmetry of cr„ and ne,i, [T^ concluded that for each a; G M, 

where F{x) is the GOE Tracy- Widom distribution function from random 
matrix theory ^0] defined in ^ below. We also find a similar result for the 
Poissonized version: 

(2) lim P I ~l <x\ = lim P | ~^ < a;! = F{x). 

We note that the length ^n(c) of the longest increasing subsequence of 
a E Inv„ has a different distribution from in. For example, while in{<^) is 
always an even integer, (.n{(^) can be both even or odd integers. Moreover, 
it was shown in [9] that 2-i(2n)^ converges to a random variable whose 
distribution function is different from F; it is given by the so-called GSE 
Tracy- Widom distribution. Hence the joint distribution of cr„ and ne^ can- 
not be the joint distribution oi in/2 and in/ 2. 
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Figure 2. The permutation matrix of the permutation a corre- 
sponding to the matching in Figure^ Since the matrix is symmet- 
ric, only the lower triangular part is shown and the entries with el- 
ement 1 are marked by x . On the left: The maximal up/right path 
(of length 2) corresponding to ne„((p((T)). On the right: The max- 
imizing down/right path (of length A) corresponding to cr„((/3(f7)) 
is realized by £g . Note that the longer down/right path indicated by 
the dashed line is not allowed as it does not fit inside the rectangles 
bounding the paths ig for any k = 1, ... ,12. 



A geometric meaning of ciCn{f{o-)) and ne„((/7(cj)) is the following. Rep- 
resent o" as a permutation matrix. Geometrically we imagine the square of 
size 2n with (1, 1) entry at the top left corner. The condition that a con- 
sists of only 2-cycles implies that the matrix is symmetric and the diagonal 
entries are zeros. Then it is easy to see that ne„((^((T)) = £„((t)/2 equals the 
length of the longest up/right chain consisting of I's in the lower-triangle 
: 1 < j < ^ < 2n}. On the other hand, for each k = 1, ■ ■ ■ , 2n, let 
^^(a) denotes the length of the longest down/right chain consisting of I's 
in the rectangle with two opposite corners (2n, 1), {k,k). Then cinificr)) 
equals the maximum of ^^(u) over k = 1, - ■ ■ , 2n [32]. See Figure l2j 



1.1. Joint distribution. The first main result of this paper is the following 
result for the joint distribution of CRt and NEj. Let F{x) denote the GOE 
Tracy- Widom distribution defined by 



(3) F{x) := exp 



(n(s) — q{s)) ds 



u{x) 



q{sfds, 
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[2n] 




Cov(cr„,ne„) 


Cor(cr„,ne„) 


4 


3 


-1/9 


-1/2 


6 


15 


-0.137777777 


-0.418918919 


8 


105 


-0.129614512 


-0.362983698 


10 


945 


-0.132998516 


-0.331342276 


12 


10395 


-0.143259767 


-0.309871555 


14 


135135 


-0.151180948 


-0.293696032 



Table 1. The exact correlation and covariance o/cr„ and ne„ for 
complete matchings of [2ri] for the first few nontrivial n 's. Note 
that both statistics are strictly negative. 



where q{s) is the unique solution of Painleve II, q"{s) = sq{s) + (?(s)^, such 
that q{s) ~ Ai(s) as s — t- oo (where Ai denotes the Airy function). The 
solution q{s) is called the Hastings-McLeod solution [32] (see also j27j). 

Theorem 1.1. Set 

(4) CR, := NE, = 

We have 

IP{CRt < x,NEt < x'} 



(5) 



' {CRt <x}f {nE, < x'} + ^'^''l^/l^'''^ + O (t-i) . 



This, together with a tail estimate, implies the asymptotics of the covari- 
ance. 

Corollary 1.1. The covariance of CUt and NEj satisfies 

(6) Gov (CRt, NEt) = \ + (t-^/^) . 
Hence, the correlation is asymptotically 

(7) p{CRt,NEt) = ^^ + 0{t-') 

where a"^ = 1.6077810345 ...is the variance of F{x) (c.f. p. 862 of [12\). 

We can also interpret CRf and NE^ as 'height' and 'depth' of certain 
non-intersecting random walks. See Section [2] below. 

We may apply the de-Poissonization argument [33J to ([s]) to find a result 
for the joint distribution of cr„ and nen,. However, intuitively, for fixed n 
and M E A^n, any G M that is used to form the maximal crossing of 

M can not be used for the maximal nesting of M. This indicates a negative 
correlation of cr„ and ne^ for a fixed n, contrary to the positive correlation 
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of CRj and NEt found in the above Corollary. This is verified for small n by 
direct computation: Table [T] shows exact calculation of the covariance and 
correlation of cr„ and ne^ for small values of n. For large n, a sampling of 
5000 pseudo-random matchings of [5000] yielded the sample covariance of 
cr25oo and iie25oo equal to —0.0420258 • • • . Therefore, a naive substitution of 
t by \/2n in ([sj) only yields the following weaker result. A further analysis is 
needed to obtain the correction terms in the asymptotic behavior of cr„ and 
ncn. A heuristic explanation for the positive correlation of the Poissonized 
random matchings is that when CRj is large, it most likely due to fact 
that the size of the matching is large and hence the maximal nesting of the 
matching is also likely to be large. 

Corollary 1.2. Set 

.o^ - _ cr„ -y/2n , _ ncn -\/2n 

^ ^ "^'^ 2-i(2n)V6' - 2-i(2n)V6- 

For each x, x' E M, 



(9) P {cr„ < X, nen <x}=F {cl-„ <x}F {n~e„ <x'} + 



\/logn 



n 



1/6 



We compare Theorem 1 1 . 1 1 with the result of [llj on the joint distribution 
of the extreme eigenvalues of Gaussian unitary ensemble (GUE). Let Amax 

(n) 

and Ajjjjj^ denote the largest and the smallest eigenvalues of n x n GUE. 
Setting 

(10) X^l := 2V2nV6(AW^ _ ^2^), x<^^ ■= 2V2,V6(^W + 
it was shown in [TT] that 

F{x'£i<xrxi:i<x'} 

= IP { aLIx < x} P {aLI < x'} + ""^^^^^"^f^f^^"^^ + O 
where Fque is the GUE Tracy- Widom distribution function defined by 



(11) 



(12) Fgue{x) := exp 



u{s) ds 



It is interesting to study the joint distribution of the extreme eigenvalues 
of Gaussian orthogonal ensemble (GOE) and compare the result with ([s]). 
This will be done in a separate paper. It might also be interesting to see if 
the error term of ([s]) can be improved to O (i~^/^) as in ( 11 ) but we do not 
pursue this in this paper. 
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1.2. Marginal distribution. We also evaluate the second order term in 
the asymptotics expansion of the marginal distributions of CR^ and NEt 
explicitly. Let [a] denote the largest integer less than or equal to a. 



Theorem 1.2. For x G M and t > 0, define xt by 

2-'^xt^/^] - t 1 

2-ltV3 + t^/3 ■ 



(13) x,:. "^r"r'-^ ' 



For each x S R, 

CRt <x} =p|NEt <x^ 
= F{xt) ^ 



(14) 



4F"(x) + \x^F'{x) 



20t2/3 

Note that since P{CRt < x} has the same value for x G [i,£ + 1) for a 



given integer £, it is natural that the leading term F{xt) of ( 194) is expressed 
in terms of Xj, which contains \t + 2~^xt^^'^\. 

In addition to this integral part correction, there is an additional shift by 
^-1/3 from X in the definition of . This is responsible for the absence of the 



term of order t in the expansion ( 194 ). For classical ensembles in random 



matrix theory, there are several papers that showed that a fine scaling can 
remove such a term (which looks like a natural term to be present.) See |24j 
for the Laguerre unitary ensemble, [37] for Jacobi unitary and orthogonal 
ensembles, [43j for the Laguerre orthogonal ensemble, and |38j for Gaussian 
unitary and orthogonal ensembles. A similar result was obtained recently for 
random growth models and intersecting particle systems in [25], including 
the height of the so-called PNG model with flat initial condition. It is well- 
known that this is precisely the length of the longest decreasing subsequence 
of random fixed-point-free involution and hence NE^. The result of [25] in 
the context of this paper is that pjNEt < x| = F{xt) + 0{t-'^l^). The 

above result finds the term of order 0{t~'^^^) explicitly. 

As in the joint distribution, the evaluation of the second order term of 
IP {crn < x} does not immediately follow from the de-Poissonization argu- 
ment in [33] . It remains an open problem to evaluate the the error terms of 
IP {cr„ < x} asymptotically. 

1.3. Toeplitz minus Hankel with a discrete symbol. Set 

oo 

(15) Gk,j{t):=Y,9k,M) 



^2n 



(2n)! 

n=0 ^ ' 
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where gk,j{n) := # {M G Mn ■ cr„(M) < A;,ne„(M) < j} so that 

oo 

P{CRt < A;,NEi < j} = Vp{crAr < k,nex < j\Af = n}F{Af = n} 

(16) 

= e-'''^Gu,,{t). 

An expHcit determinantal formula of Gk,j{t) was obtained in [15] which we 
describe now. 

Stanley had shown earlier that matchings are in bijection with oscillating 
tableaux of empty shape and of length 2n (see section 5 of [E]). This 
was further generalized to a bijection between partitions of a set and so- 
called vacillating tableaux in jl5j . In the same paper, it was shown that the 
maximal crossing (resp. nesting) of a partition equals the maximal number 
of rows (resp. columns) in any partitions appearing in the corresponding 
vascillating tableau. 

Since an oscillating tableau can be thought of as a walk in the chamber 
of the affine Weyl group C„, gk,j{n) equals the number of walks with n 
steps from — 1. . . . , 2, 1) to itself in the chamber < Xj < • • • < X2 < 
xi < j + k + 1 where each step is a unit coordinate vector or its negative in 
V . The number of such walks was evaluated by Grabnier in [31] using the 
Gessel-Viennet method of evaluation of non- intersecting paths. This result 
implies (see the displayed equation before (5.3) in [T5]) that 

" 1 2m-l / , \ "1 

(17) Gk,j{t) = det 



1 sin f^) sinful e2*™^(WH 
- — =0 ^ ^ 



a,b=l 



where 

(18) m:= i + k + l. 

We prove Theorem |1.1| by analyzing the determinant (17) asymptotically. 
For this purpose, we first re- formulate the determinant slightly. By writing 
the product of the sine functions in terms of a sum of two cosine functions 
and noting the realness of the entries, we find that 

(19) {t) = det [ha-b - /ia+fe]if,=i 
where 

2m- 1 



(20) hi := — ^ g-i7rrf/mg2tcos(7rr/m)_ 



2m 

r=0 



This is the determinant of a Toeplitz matrix minus a Hankel matrix. This 
structure is important in the asymptotic analysis. An interesting feature of 
the above determinant is that the measure for the Toeplitz determinant is 
not an absolutely continuous measure but a discrete measure. 
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Let u) := e'^*/™ be the primitive 2m*'' root of unity. Define the discrete 
measure 

2m- 1 

(21) d^x^{z) ■=,^Y. e*^'+'"'^<5-'-(^) 



2m 

r=0 



on the circle. Let 'Kn,m{z) be the monic orthogonal polynomial of degree n 
with respect to d^m-, defined by the conditions 



(22) i) Z 'Knm{z)dflm{z) = 0, < £ < n. 

J\z\=l 

We emphasize the dependence on m since later we will use the notation vr^^oo 
to denote the case when 'm = oo'; the orthogonal polynomials with respect 
to the absolutely continuous measure e**^^"*"^ ^2^' ^o^e that dfj^m depends 
on the parameter t and hence 7Tn,m (z) also depends on t. When we wish to 
emphasize this dependence on t, we write ■Kn,m{z]t). 

The fact that the t-dependence of the measure is from the factor e**-^^^ 
implies the following basic formula, which is proved in Section [3] below. 



Recall from ( 18 ) that m := j + k + 1. 
Proposition 1.1. We have 
(23) logP{CRt < fc,NEt < j} 

ft f S 

= / vr2,+i,„(0; T)dT+ / / Qf{T)dTds 
Jo Jo Jo 

where 



|2\ 



(24) Qf{T) := -{TT2j,m{0; r)7r2j+2,m(0; r) + k2i+i,m(0; r) 

+ 7r2j,m(0; r)7r2j+2,m(0; r)|7r2j+i,m(0; r)p. 

We obtain the asymptotics T^2j+i,mi0, t) for ^ = 0, 1, 2 by using the asso- 
ciated discrete version of the Riemann-Hilbert problem (see e.g. [6]). See 
Sections |4j [5] and |6] below. 



We compare the analysis of this paper based on the formula (23) with 
the analysis of the determinant of a similar Toeplitz minus Hankel matrix 
in [9]. Even though the determinant in [9] was for continuous measure 
(which is precisely the one for the marginal distribution of NEf ; see Section [s] 
below), the basic structure of the matrix is the same; a Toeplitz minus a 
Hankel matrix. Denoting the matrix by Dj, the approach of [9j was to write 
Dj = Doo Y[^=j d'I^ where Doo is the strong Szego limit, which exists in 
that particular case, and analyze Dn/Dn+i, which can be evaluated from the 
Riemann-Hilbert problem for the n*'* orthogonal polynomial. For our case, 
since the measure is discrete, the strong Szego limit does not apply. Indeed 
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Dn = for all large enough n. Then alternatively one can still analyze Dj by 
expressing Dj = Dq Y{n=i TT^ ^ done in [3|. However, this expression 
is more subtle to analyze since log(L'„/-Dn4-i) is not small when n is small 
(indeed it grows as n decreases when t is proportional to j) and this requires 
careful cancellations of the terms in the product. Though this was done for 
the leading term in |^ , the evaluation of the lower terms in the asymptotic 
expansion in this method becomes more complicated. A particularly useful 



point in using formula ( 23 ) is that we only need to consider the so-called full 
band case (and the transitional case when a gap and a saturated region are 
about to open up) in the Riemann-Hilbert analysis. This makes the analysis 
much simpler and it becomes easier to evaluate the lower order terms. On 
the contrary, if we use the expression Dj = -Dq nn=i iJ^' then we need 
to consider both the so-called void-band case and the saturation-band case, 
including the transitional cases, in the Riemann-Hilbert analysis (and this 
is the reason for the need of cancelations mentioned above.) 

The continuous Riemann-Hilbert problem for 7r„^oo(-2;i) was analyzed 
asymptotically to the leading term in [H El [3]. We expand this work to 
the discrete counterpart and moreover, we improve the analysis so that 
we compute explicit formulae for the first three terms in the expansion 
of the solution in both the discrete and continuous cases. As a technical 
note, we remark that we use a different local map for the so-called Painleve 
parametrix related to the local problem for the Riemann-Hilbert problem 
from the previous cases [H [19] . We adapt the map used in the recent paper 
|14j for a different parametrix, which seems to be useful for further analy- 
sis in other Riemann-Hilbert problems. For the purpose of this paper, we 
only analyze the full band case (and the transitional case) of the discrete 
Riemann-Hilbert problem. The analysis for the full parameter set of the 
discrete Riemann-Hilbert problem will be discussed somewhere else in the 
context of Ablowitz-Ladik equations and Schur flows in integrable systems. 

A determinantal formula of the marginal distribution P {NE^ < j} can be 
obtained from the joint distribution by taking — )• cx) while keeping j fixed. 
Then we find a Toeplitz minus a Hankel determinant with symbol e*^^"*"^ . 
Here too, the factor of e**-^"*"^ in the limiting measure implies a formula 
for the marginal distribution analogous to (23). See Section [3] below. 



The Toeplitz determinant with symbol e**-^"*"^ is known to be describe 
the distribution of the length of the longest increasing subsequence of a 



random permutation |2S|. By using a formula similar to (23), the analysis 
of this paper implies the following result. 

1.4. Longest increasing subsequence of random permutation. Con- 
sider the symmetric group Sn of permutations of size n and equip it with the 
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uniform probability. Let /^(vr) denote the length of the longest increasing 
subsequence of n G Sn- Let A/i be a Poisson random variable with param- 



eter and let Lt denotes It was shown in |3| that 



-2t 



tl/3 



converges to 



the GUE Tracy- Widom distribution (12). We evaluate the next term of the 
asymptotic expansion explicitly. 

Theorem 1.3. For each x G M, 



(25) 



where 



(26) 



Lt-2t 
ti/3 



< X 



Fgue{x 



1 



it)) 



10*2/3 



Fgue{x) + ^x^F'(jue{x) 



+ o (t-i) , 



X 



it) 



[2t + xt^/^] - 2t 



The study in [53] also considered the height of the so-called PNG model 
with the droplet initial condition, which is distributed precisely as Lt, and 
showed that the above distribution function is Fgue{x^^^) + 0(t~2/3)_ ^he 
above theorem evaluates the error term explicitly. 

For the Gaussian unitary ensemble, [T71 [18] evaluated the distribution 
of the largest eigenvalue explicitly up to the term of order 0(n^2/3') -v^j^ich 
corresponds to the term of order i~2/3 ^^le above expansion. It would be 
interesting to compare the term in the above theorem with the formula of 

[r7i[iE]. 



1.5. Organization of paper. In Section [2j we consider a non-intersecting 
random process that gives rise to CKt and NE^. Proof of Proposition [3] is 
given in Section [3| The Riemann-Hilbert problem is introduced in Section [4j 
and is analyzed asymptotically in Sections [5] and [6j Theorem |1.1| and Corol- 
lary [TTT] are proved in SectonjTJ and Theorem 1 1 . 2 1 and Theorem |1.3| are proved 
in Section [8j We prove Corollary 1.2 in Section [9] using a de-Poissonization 
argument. Finally, the Riemann-Hilbert problem for the Painleve II equa- 
tions that are needed to model the local parametrix of the Riemann-Hilbert 



problem for orthogonal polynomials are discussed in Section 10 
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2. Height and depth of non-intersecting continuous-time 

simple random walks 

In Section [l.3| we discussed a relation between cr„ and ne„ and a walk in 
the chamber {0 < < • • • < X2 < xi < j -|- /c + 1} of the affine Weyl group 
Cn- In this section, we give an interpretation of CRt and NEt in terms of 
the 'height' and 'depth' of continuous-time simple random walks. 

Let N^{t) and N~{t) be two independent Poisson processes of rate 1 
and set Z{t) := N^{t) — N~(t) be a continuous-time simple random walk. 
Then Z{t) is an Z-valued Markov process with the transition probability 
Ps{a,b) = ps{a - b) where ptia) = e~^^Y.n& = a e Z. 

Here we used the convention that l/nl = if n < 0. Set 

(27) cPiz) := ^ {e-^'ptia)) z'^ = e'^^+^-'\ 
Then we have 

(28) pt{a) = e-2V_a = e^^^^, (j^a := i z^Xz) — . 

J\z\=l 

Let Zi{T), i = 0,1,2,..., be independent copies of Z{t), and consider 
the infinite system of processes Xi{T) = Zi{T) — i, i = 0, 1, 2, . . . . Fix a 
number t > 0. We will consider the process conditioned on the event that 
(a) Xi{t) = Xi{0) for all i and (b) Xi^r) do not intersect in time [0,t], 
i.e.Xo(r) > Xi{t) > ... for all r G [0,t]. A precise interpretation will be 
given below. Such non-intersecting continuous-time simple random walks 
have been studied, for example, in [lH [2l 

Define the 'height' K := max^-gjo,*] Xo(t). and define the 'depth' J as the 
smallest index such that Xi{T) = i for all r G [0, t] and for all i = J, J+1, .... 
In other words, only the top J processes moved in the interval [0,t]. We 
are interested in the joint distribution of J and K conditional of the above 
event satisfying (a) and (b). 

Precisely, fix G N and let SIat and ^Bat be the events defined as 

(29) := {Xi{t) = Xi{0) = -i, i = 0,l,...,N-l}, 

(30) 5Sjv :={Ao(r) > Ai(t) > ••• > Ajv-i(t) > -A + l, TG[0,t]}. 

The condition that A7v_i(t) > —N + 1 for all r G [0, t] is natural because 
J is likely to be a finite number and by definition of J, Xj_i(r) > Aj_i(0) 
for all r G [0,t]. The joint distribution of K and J is interpreted as 

(31) P{k,j) := lim F (K < k, J < j D ^n) ■ 

N^oo 
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Lemma 2.1. Let K and J he the 'height' and 'depth', respectively, defined 
above. Then 

(32) P(A;,i) = e-*'/2Gfc,,(t) 



where Gkj{t) is given in (19). 



Proof. We first evaluate P(2tAr n QSat). The condition that Xi(r) > -N, i = 
0, • • • , A^— 1, implies that (t) has an absorbing boundary at — A^. Since the 
transition probability of Xi with an absorbing boundary at — is pt{a, b) — 
Pt{—2N — a,b), the Karlin-McGregor formula [40j of non-intersecting proba- 
bility applied to continuous-time simple random walks (see e.g.[Il[2j) implies 
then that 

P (2tiv n 5Sjv) = det [pt(-a, -6) - pt(-2iV + a, -6)]^fe-io 
= e-^*^det[0._,-0.+,]^,^,. 

Secondly, we evaluate P {{K < k, J < j} H SIat n ^n)- We assume that 
is large so that N > j. By the definition of K and J, the desired 
probability equals P(Cn !l)) where C and Ti are independent events defined 
as follows. C is the event that the top j processes, Xo(t), . . . , Xj_i(r), 
satisfy the two conditions (a) Xi{t) = Xi{0) for all i = 0, . . . , j — 1 and 
(b) -J + 1 < Xj_i{T) < ••• < Xo(t) < k for all r G [0,t] i.e.the j non- 
intersecting paths are not absorbed at the boundaries —j and k + 1. T) 
is the event that Xi{T) = —i for all i = j, j + 1, . . . , N — 1 and for all 
r G [0,t] i.e.the bottom N — j processes stay put during the interval [0,t]. 
Clearly, P(2)) = (e~^*) "'. On the other hand, from the Karlin-McGregor 
formula again, P(C) = det [pt(— a, — 6)];^'^^g where pt{a,b) is the transition 
probability of Z{t) in the presence of the absorbing walls at —j and k + 1 
in time t. It is easy to see that 

(34) pt{a, 6) = ^ [pt{a + 2nm, b) - pti-2j - a + 2nm, b)] 



where m := j + k + 1. Now consider the identity z~"(/)(z) = X^nez 

Set CO := e'^*/"*. By inserting z = to'' , r = 0,1, ■ ■ ■ , 2m — 1 and summing over 

r, we find that 

2m- 1 

(35) i^T^H^l = 2mJ2 4>a+2mn, ^ := e-^/'". 

r=0 neZ 



Hence from (28), (34) becomes 

(36) ptia, b) = e~^* {ha-b - h_a-b+2j) 
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where 
(37) 



ha 



\z\ = l 



^ 2m-l 



2m 



r=0 



Hence, for N > j, 

(38) F{{K <k,J <j}n^Nr]^N) =e~ aei[ria-b 

The strong Szego hmit theorem for Toephtz minus Hankel determinants 



-2*^det [K-b-ha+b]ib=i- 



(see for example |10j ) imphes that for the function (^{z) in (27), det 



'}a-b 



i>a+b]a,b=l 



e* as — )• oo. Therefore, from (33) and (pSl) we find that 



(39) P{j, k) = hm 



det [ha-b - ha+b]i 6=1 



N^oo det [(f>a-b 



k 1^ 



e det [ha-b - ha+bll 



b=l ■ 



This is (32). 



□ 



Hence K and J have the same joint distribution as CR* and NEj. This 
non-intersecting process interpretation of CRj and NE^ provides some useful 
information. As an example, note that the process considered above has a 
natural dual process (see Figure pi). In the dual process the roles of K and 



Figure 3. A non-intersecting continuous-time simple random 
walks (left) and its dual walk (right). 

J are reversed: the depth is K and height is J in the dual process. It follows 
that K and J, and hence CRj and NE^, are symmetrically distributed. 

In various non- intersecting processes, including the above model, the top 
curve is shown to converge, after appropriate scaling, to the Airy process 
in the long-time, many- walker limit (see e.g. |34| I36j). Then it is natural 
to think that the leading fluctuation term of K is given by the maximum 
of the Airy process. It is a well-known fact that the maximum of the Airy 
process is distributed as the GOE Tracy- Widom distribution. This was first 
proved indirectly in [SS]. A direct proof was only recently obtained in [2D]. 
(See also [H] for the distribution of the location of the maxima.) Hence the 



leading term F(x) in (194) is as expected. Moreover, when t becomes large, 
it is plausible to expect that the fluctuation of the top curve of the original 
process (whose max is K) and the fluctuation of the bottom curve of the dual 
process (whose min is — J) become independent at least to the leading order. 



14 



JINHO BAIK AND ROBERT JENKINS 



The leading term of Theorem 1 1 . 1 1 is natural from this. Theorem 1.1 evaluates 



the second term of the asymptotic expansion of their joint distribution. 

For a family of finitely many non-intersecting walks, it is interesting to 
consider the maximum of the top curve and the minimum of the bottom 
curve. It is curious to check if the joint distribution of them would have 
the same expansion as in Theorem |1.1[ This will be considered elsewhere. 
Finally, we mention that the asymptotics of the distribution of the width of 
non-intersecting processes was studied recently in [7]. 



3. Proof of Proposition [TTT 



In this section, we give a proof of Proposition |1.1[ We also obtain sim- 
ilar formulas for the marginal distributions of CR^ and NE^, and for the 
distribution of Lf. They are stated at the end of this section. 

Let dp be a (either continuous or discrete) measure on the unit circle and 
define a new measure dp{z; t) which depends on a parameter t as 

(40) dp{z;t):=e'^'+'~'^dp{z). 



The measure (21 ) associated to the joint distribution of CRt, NEt is certainly 
of this form, but the following algebraic steps apply to general dp. 
Let 



(41) he{t) := 6 z-'^dp{z;t). 

J\z\=l 

We are interested in finding a simple formula for the second derivative of 
the Toeplitz determinant r„(t) and the Toeplitz - Hankel determinant Hn{t) 
(see ([Tol)) associated to the measure dp{z;t): 

(42) Tn{t):=det[ha-b{t)]l,^,, F„(t) = det [/ia_fe(t) - /ia+,(t)]",^^. 

We assume that when dp is a discrete measure, n is smaller than the number 
of points in the support of dp. 

Let TTniz; t) = z"" + ■ ■ ■ , n = 0, 1, 2, • • • , be the monic orthogonal polyno- 
mials defined by the conditions 



(43) < TTn, z'- >:= 6 7r„(z; t)z'^dp{z] t) = 0, < i < n. 

J\z\ = l 

Set 

(44) Nn{t) : = < TTn, 7r„ > = < TTn, > . 

Then it is well-known that (see e.g. section 2 and 3 of [8j for the second 
identity) 

(45) T,{t) = H Nr^it), H,{t) = II iV2n(i)(l - vr2„(0;t))-i. 

n=0 n=l 
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Define (see !39l) 



(46) <(z; t) := z"7r„(z-i; t) = 1 + a^z + ■■■+ alz''-^ + 7r„(0; 
This polynomial satisfies the orthogonality properties 

(47) (^ttI, z^^ = NnSkfi, A; = 0,l,...,n. 
Recall the Szego recurrence relations |49j : 

■Kn+l{z) = ZTTn{z) + 7r„+i (0)7r* (z) , 

ZTTn[z) = — (vr„+i(z) - 7r„+i(0)7r„+i(2;)j . 

The second relation, when we compare the coefficients of z^~^^, gives rise to 
the relation 

(49) %^ = l-K+i(0)|^ 

We now derive a differential equations for 7r„(0;t) and Nn{t). All the 
differentiations are with respect to t and we use the notation /' for ^f. By 
differentiating the formula (tt^, z^^ = 0, k = 0, . . . ,n — 1, we obtain, by 
noting ^e*(^+^~') = (z + z-i)e*(^+^"'), that «, z'=) + (7r„, z''+^ + z^-'^) = 
0. Then by using the orthogonality conditions, we find that 

(<, z'=)=0, fc = l,...,n-2, 

(50) «, 1> = - (7r„, Z-^) = - {zTTn, 1) = ^n+l(0)iVn, 

where the last equality in the second condition above follows from the first 
recurrence in (|48|). From these relations, we conclude that, for n > 1, 



(51) 7r'Jz;t) = -M) (,r„+i(0;t)<_i(z;t) -7r,_i(z;t)). 

This can be checked by taking the difference and noting that the difference is 
a polynomial of degree at most n— 1 and is orthogonal to z^ , k = 0,1, ■ ■ ■ ,n— 



1. Evaluating ( |5l| at z = 0, we obtain, using (49), for n > 1. 

(52) 7r;(0; t) = (^„+i(0; t) - 7r„_i(0; t))(l - |7r„(0; t)\^). 

This equation is related to the Ablowitz-Ladik equations and the Schur flows 
(see, for example, [35ll30] ). 

We also differentiate Nn{t) = {-Kn, vr„) and obtain 

(53) K = 2 «, vr„) + 2 (z7r„, 7r„) =< 2z7r„, 7r„ > . 



Using the first recurrence of ( 48 ) , 

(54) {zTTn, TTn) = (vr„+l, 7r„) -7r„+l(0) {-Kn, VT*) = -7r„+i (0)7rn (0) {iTn, z"^) . 
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Hence, we obtain, for n > 0, 

(55) Kit) = -2^„+i(0; t)TTn(,0; t)Nn{t). 



We now evaluate the logarithmic derivatives of Tj and Hj. From (45) 



and (55), we find that 
(56) 



(iogr,(t))' = ^ 



n=0 



Nnjt) 
Nn{t) 



i-1 

-2^^„(0;t)7r„+i(0;t). 

n=0 



We take one more derivative. By using (52), for n > 1, 
(57) (7r„(0)7r„+i(0))' = P„+i-P„, 

where P„ := |7r„(0)|2+7r„_i(0)^„+i(0)(l-K„(0)|2). Forn = 0, (7ro(0)7ri(0))' = 
7r'^(0) = (7r2(0) — 1)(1 — |vri(0)p) = Pi — 1. Hence from a telescoping sum, 
we obtain 

(58) 

^ (log (e-*V,(t)))" = -(7r,_i(0)7r,+i(0) + K,(0)|2)+vr,_i(0)7r,+i(0)|7r,(0)p 



We now consider Hj{t) in (45). By taking the log derivative and using 



(52), (55) and 7ro(z) = 1, 



(59) 



{\ogH,{t))' = Y, 



n=l 



N2n 1 - vr2„(0) 



2i 



vr2i+i(0)- J];7r„(0)7r„+i(0). 



n=0 



From (56), we find that 
(60) 



{\ogH,{t))' = 7r2,+i(0) + - (logr2,+i(t))' . 



Proposition is proven from (16), (19), (58), and (60) by noting that 
7r„(0; 0) = for all n > 1, and T, (0) = 1 and ^^, (0) = 1 for all j > 1. 

The marginal distribution of NE^ is obtained from (|16|) by taking the 



limit k — )• oo. Then by taking m — )• oo in (19), we find that P{NE( < j} 



e * ^'^Gooj where Goo,j{t) is same as (19) where the measure in (21) is 
replaced by 

(61) <i,^(„:.,*.-)_^!_. 



Then the above computation applies that 

(62) logP{NEt< i}= / ^2,+i,oo(0; r)dr+ / / Qf(r)drds, 

Jo Jo Jo 
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where TTn,ooiz;t)) is the monic orthogonal polynomial of degree n with re- 
spect to the measure (61 ) and Q°°(r) is same as (43 ) with 7rn,m{z; r) replaced 
by vr„,oo(2; r). Due to the symmetry, P{CRt < j} = P{NEt < j}. 

Finally, it is well-known, |29| H6] , that for the length Lt of the Poissonized 
random permutation defined in Section 1.4, ¥{Lt < i} = e~* Ti{t), where 



Tj (t) is the determinant of the £x£ Toeplitz matrix ( 42 ) with respect to the 
meas 

(63) 



measure (61). Hence we have 

log F{Lt<£} = 2 f r QfU (r) dr ds. 
Jo Jo 2 



4. Orthogonal Polynomial Riemann-Hilbert Problems 

We prove Theorems |1.1| and 1.2 by deriving asymptotic expansions of 
vrn,m(0; r) and 7rn,oo(0; r), n = 2j, 2j -|- 1, 2j -|- 2, r G (0, t), in the joint limit 
t, j, m — )• cxD such that given any fixed x, x' G M: 



(64) 



J 



t + -t 
2 



2; -1/3 



k = t^ 1 



X .1/3 



m 



j + k + l. 



The jumping off point for our analysis is the fact that iTn,m{z; t) and nn,oo{z; t) 
can be recovered from the solution of the following discrete and continuous 
measure Riemann-Hilbert problems respectively. 

Riemann-Hilbert Problem 4.1 for discrete OPs Find a 2 x 2 matrix 
'y{z; t, n, m) with the following properties: 

1. Y{z;t,n,m) is an analytic function of z for z G C\{uJr}'^^Q^ where uj^- := 
w'' and oj := e*'^/"^. 

2. Y(z; t, n,m) = [I + (l/z)] z"'^^ as z ^ 00. 

3. At each ujr, Y{z;t,n,m) has a simple pole satisfying the residue relation 
(65) Resz=(^^ Y(2;; t, n, m) = lim Y(z;t, n) 



2m 



As is well-known (see e.g. |27j . [6j), and may be verified directly, the 
solution Y(z; t; n, m) is given by 

■^n,m(^) t) * 
-<-l,m(^;0/^n-l,m 



(66) 



'Y{z; t; n, m) 



where we recall that tt* ,^ is the reverse polynomial defined by (46) and 

,(a;,;t)w-"+ie*(^''+^"') 



Yi2{z;t,n,m) 
Y22iz;t,n,m) 



^ 2m- 1 
1 sr^ TTr, 



2m 



r=0 



Z — IjJr 



2m- 1 



2m ^ 



r=0 



Z — UJr 
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Hence, using the OP properties listed in (43)- (48) we can easily check that 
(67) Y(0; t,n,m) = 



"17,771 



7r„,m(0) N„ 

(0). 

Note that the generic (2, 2)-entry would be 7r„,m(0) but as our weight e*^^"*"^ 
is real 7r„,m(0) = vr„,m(0). 

The continuous RHP can be thought of as a limit of the discrete case when 
m, the number of points in the support of the measure, goes to infinity. 

Riemann-Hilbert Problem 4.2 for continuous OPs Find a 2x2 ma- 
trix Y°°(z; n) with the following properties: 

1. Y°°(z;t, n) is an analytic function of z for z € C\S, Ti := {z : \z\ = 1} 
oriented counterclockwise. 

2. Y°^(z;t,n) = [I + O {I / z)] z"""^ as z ^ oo. 

3. Y°° takes continuous boundary values YJ^ and YJ^ as z — )• S from the 
left/right respectively satisfying the relation 

Y~(z;t,n)=Y-(z;t,n) (q % j , z G S. 

The solution Y°° is related to the orthogonal polynomials vrn^oo with re- 



spect to the measure /ioo (61) and we have 

7r„,oo(0) 



Y'^iO;t,n,m, 

\-L/I\n-l,oo 7r„^oo(Uj 

Precisely this continuous Riemann-Hilbert problem was analyzed asymp- 
totically in [H El [6] . The steepest-descent analysis for discrete Riemann- 
Hilbert problem was studied for general discrete measure on the real line in 
[6] . Both works expand upon the continuous weight case studied in [211 ES] • 
In the course of proving Theorems |1.1| and 1.2 we improve these results as 



follows: we expand the analysis of [6j to the case when a gap and saturated 
region of the equilibrium measure (see the discussion below) are about to 
open up, and we compute explicit formulae for the first three terms in the 
expansion of the solution in both the discrete and continuous cases extending 
the results of [H |8l [6] where only leading terms were calculated. 

One of the key steps in the steepest-descent analysis of Riemann-Hilbert 
problems is the introduction of the so-called (^-function. For the Riemann- 
Hilbert problem |4.1| for discrete orthogonal polynomials, the (^-function is 
given by g{z) = J^^^^i log{z— s)dfj,{s) where dfi{s) is the so-called equilibrium 

measure satisfying < d^{s) < (^^^ ^-S- El)- '^^^ upper-constraint 

diJ,{s) < *° ^^^^ weight is discrete. The support of 

d^ consists of three types of intervals, voids (where dfj, = 0), bands (where 
< dfiis) <^^J and saturations (where d/x(s) = ^^,). 
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For the continuous Riemann-Hilbert problem, the upper-constraint for the 
equihbrium is not present and there are no saturations. For the Riemann- 
Hilbert problem 4.2, it was shown in [4J that with 7 = the equilibrium 



measure is supported on the entire unit circle, i.e., S is a single band when 
7 > 1, and S consist of single void and band intervals with the void set 
centered about z = —1 when 7 < 1, . 



In the discrete problem 4.1 the solution Y now depends on the three pa- 
rameters (t, n, m) and as we shall see in Sectionjsj the equilibrium measure's 
support depends critically on the two parameters 

n 2m — n 

(70) 7 = — and 7 = . 

^ ' ' 2t ' 2t 

As each of these parameters passes through the critical value ^crit = 1 a 
transition occurs in the support of the equilibrium measure. 

It turns out that to prove Theorems 1.1- 1.3| we only need to evaluate 



Y(0; t,n,m) in two regimes: the 'exponentially small regime' 

(71) n>2t{l + 5), 2m - n > 2t{l + 6) 
for a fixed 6 > 0, and the 'Painleve regime' 

(72) 2t- Lt^/^ <n<2t{l + 6), 2t - Lt^'^ <2m - n <2t{l + 5) 

for fixed L > and (5 > 0. In the 'exponential' case 7, 7 > 1+5 and the equi- 
librium measure is supported on the whole of S, while in the 'Painleve' case 
7,7 G [1 — 1 + (5] and the equilibrium measure is in the transitional 

region where a void and saturation region are beginning to open at z = — 1 
and z = 1 respectively. As such we never need to consider cases in which 
either a void or saturation have fully opened and we restrict our attention 
to the full band (and the transitional) case only, focusing on obtaining the 
three lower-order terms of the asymptotic expansion explicitly. In this case 
the (7-function is explicit and the transformations of the Riemenn-Hilbert 
problem will be all stated explicitly without mentioning the (7-function in 
the subsequent sections. 

There are many interesting related problems in which one needs an as- 
ymptotic description of the iTn^m for a whole range of degrees n; one such 
example which we plan to study in the future is the Ablowtiz-Ladik equa- 
tions. There we will fully describe the structure of the equilibrium measure 
in the full range of parameter space. 

The analysis of the discrete and continuous Riemann-Hilbert problems 
have strong similarities and we analyze them simultaneously. The important 



^This is actually the inverse of the parameter appearing in 4 which we find more 
convenient to work with presently. 
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fact, which we clarify in Sections 6.2 6.3 is that in the discrete Riemann- 



Hilbert problem we can partition the solution into terms that come from 
(two) continuous Riemann-Hilbert problem which correspond to the mar- 
ginal distributions and the remaining 'joint' terms which contribute only to 
the joint distribution. 

5. The Exponentially Small Regime 

The first steps of the steepest-descent analysis are the same for both 
the exponentially small regime and the Painleve regime. We begin by first 



considering parameters {n,m,t) in the 'exponentially small regime' (71): 

n>2t{l + 6), 2m-n>2t{l + 6) 

for fixed 6 > 0. We assume that 6 < 1/2 (See the discussion before ([90|)). 

We begin our analysis of RHP |4.1| by first introducing a transformation 
Y I— )• Q such that the new unknown Q has no poles. Let E denote the unit 
circle and let Si„ and T,out denote positively oriented simple closed contours 
enclosing the origin such that Sj„ C {z : |2:| < 1} and ^out C {z : \z\ > 1}; 
let and Q- denote the non-empty open sets enclosed between S and Sj„ 
and S and T,out respectively. Define 




Figure 4. The contours and regions used to define the map Y i— > 
Q. The contours and Sou* can be deformed as necessary pro- 
vided they do not intersect S. 



(73) 



Q(^) 



Y(z) 



Yiz) 



1 



z G 17- 



The triangular factors introduced in the above definition have poles at each 
(jOr and the residues are such that the new unknown Q(z) has no poles, 
but is now piecewise holomorphic. Note that the residue of each triangular 
factors at each z = ujr are the same since z^*" = 1 at z = cjr. Two different 
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extensions of Q as above were introduced in [39] (see also 0). By explicit 
computation Q,{z) satisfies 

Riemann-Hilbert Problem 5.1 for Q(z): Find a 2 x 2 matrix Q(z) such 
that 

1. Q(z) is analytic in C\(S U Tim U Tout)- 

2. Q(z) = [I + (1/z)] z"'"^ as z^oo. 

3. Along each jump contour Q+(z) = Q_(z)Vq(z) where 



' ( 1 _2~"e*^^'''^ '^^^ 



(74) 



1 



,0 f 



rZ e 



1 1 ^-ra„t(2:+2 ^) 

-■^ ^2™ — 1 



z G S 



z € S 



out 



Once we transforms a RHP with poles to a 'continuous' RHP as Q, the 
next step is to introduce a "(7- function" . However, for the above RHP, when 



the parameters are in the regimes (71) and (72), it turns out that the g- 



function is simple and explicit. We proceed by explicitly defining 



(75) 



S(z) := 




\z\ < 1 



\z\ > 1 



Clearly Y(0) = S(0) and S(z) = I + O (z'^) for large z. Calcu- 

lating the new jump matrices, we arrive at the following problem for S(z). 

Riemann-Hilbert Problem 5.2 for S{z): Find a 2 x 2 matrix-valued 
function S(z) such that 

1. S(z) is analytic for z £ C\(S U Sj„ U Tout)- 

2. S(z) = / + O (1/z) as z^ 00. 
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3. The boundary values ofS{z) satisfy the jump relation S+{z) = S-{z)Vs{z) 
where 



(76) Vsiz) = { 




zgT. 



z 



zGT, 



out 



where 



(77) 



0{z-n) :-- 



{z-z )+7log(-z), 7: 



{z- z ^) - 7 log z, 



7 



n 

2m — n 
2i ■ 



Here the log is defined on the principal branch. 

Now we assume that the parameters are in regime (71). Note that for 
G S, 0(e*°) G iM. Also note that writing z = re*", we have 
^ = cosa+7 > -1+7 > 5 > Oand 1^ [Re6i(re*";7)]^^-^ = 
7r~^ < — 7r^^ < if r > 7^^. Hence Re0(re*";7) < 
(— l+7)(r — 1) for r G (7""*^, 1) and for all a G (— vr, vr]. Therefore, for a given 



any e 



^[Re%e-;7)], 



-r^ cos a 



6 > 0, there exist < ri < r2 < 1 and c > such that Re 



leire'";j) 



< 



for all r G [ri,r2], a G (— tt, tt] and for the parameters {n,m,t) in the 



regime (71). Note that this implies that 

nRe[ie(re'";7)] 



(78) 



=2t9(^;7)| 



< e" 



for parameters (n,m,t) in the regime (71) 



Similarly, Re ^0(-e*°;7) < — c for all r G [?'i,r2], a G (— 7r,7r] and for 



the parameters (n, m, t) in the regime ( 71 ). This can be easily seen by noting 
that (l){z;j) = 6{—z~^;^). Hence 



(79) 



= 2t0(2:;7)| 



g(2m-n)Ro[ie(re'";7)] ^ g-c(2m-n) — < l^l < — 



ri 



^2 



for parameters {n,m,t) in the regime (71). 

Let Cin-i, Cin^i, Cout,!-, and Cout,-i be the contours as depicted in Figure 
[5] such that Cin-i and Cin,i lie in the annulus ri < \z\ < r2 and Cout.i and 
Cout-i he in the annulus ^ < |2;| < ^. Make now the following change 
of variables which moves the oscillations on S into regions of exponential 
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decay. 



(80) T{z) = { 



S{z 

S{z 

S{z 
[S{z 



\ (-l)"e2*«\ 
.0 1 ) 

.0 1 ) 

' 1 0^ 

^(-l)"e-2*^ \, 

' 1 0^ 



1 0^ 

.g-2t<^ 1 



.0 



z G 

z G 0_,o 

z e 

elsewhere 



Note that Y(0) = T(0) |^ ^ . ExpUcitly calculating the new jumps, the 
new unknown T{z) satisfies the following problem 

Riemann-Hilbert Problem 5.3 for T(z): Find a 2 x 2 matrix-valued 
function T{z) satisfying the following properties 



1. T{z) is analytic in C\(Si„ U T^out U Cm,±i U Co„t,±i). 

2. T(z) = I + 0{l/z) as z^oo. 
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3. The boundary values ofT{z) satisfy the jump relation T^{z) 

where 
(81) 

Vt{z) = 



T.iz)VTiz) 



1 -(-l)"e^*'' 

1 

1 
(-l)"e-2*» 1 

1 

1 e^*^ 
1 



Z G Cin, — 1 

Z ^ Cout,-l 

Z ^ Ciri,l 

Z € Cout.l 





1 0^ 






1- 


ly 




-(-!)"<= 





1 







Z ^ Sin^ — l 



(1 - z2m.)-<T3 ^ g J., 



(-1)"'= 



1 0^ 

2te I Z (z Sou£,l 



1 ^2i0 



Then from ([TsJ) and ((79]), we find that Vt{z) = I + O (g-^'^^^'t'^'^'"-"}) 
uniformly for z on the contour. Hence we obtain the following result. 



Proposition 5.1. Let Y{z;t,n,m) be the solution to the RHP (4.1). For 

any 6 > 0, there exists a constant c > such that, if 



(82) 
then 



(83) 

In particular, 
(84) 



n>2t{l + 6), 2m-n>2t{l + 6), 
Y(0; t, n, m) (^^ = I + O (^g-^^^^^i^'^"-"}) . 



vrn,m(0;t) = O (^e-=max{n,2m-n}^ 



6. Painleve Regime 



We now consider the parameters {n,m,t) in the regime (72): 



(85) 2t- Lt^^^ <n<2t{l + 6), 2t - Lt^'^ < 2m - n < 2t{l + 5) 

for fixed L > and 5 > 0. We assume that 5 < 1/2 (see the discussion 
before ([90j)). 

Let S(z) be same as in the previous section. When 7 e [1 — 5, 1 + 5], 



the estimate (78) does not hold any more. However, it is easy to check 
using a similar calculation as before that the exponential decay still holds 
in an annular sector away from the point z = —1. (Note that the (double) 
critical point of 0{z] 1) is z = —1.) More precisely, one can check that given 
5 £ (0, 1), there exist positive constants ai > 0((5^/^) > and pi > 0(6^/'^) 
such that if 7 G [1 — 5, 1 + 5], then |e^^^^''''''| < 1 for z in the annular sector 
Sin-i ■= {z = re*" : pi < r < 1, |a| < vr — ai}. Moreover, if z is in a 
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7=0.99 




7=1.01 




Figure 6. The sign of I{.e6{z;j) for values of ^ near jcrit = 1- 
Note the sign change near z = ~1 on either side of the transition. 



compact subset of 5m, -i, then there exists c > such that [e^^^^'"^)! < 
uniformly in 7 G [1 — (5, 1 + 5], see Figure |6j 

Similarly, from the symmetry 4>{z]'j) = 6{—z~^;'y), under the same as- 
sumptions, |e^''^*^^'^^| < 1 for z in the annular sector Sout,i '■= = ^e*" : 
pi < r < 1, ai < |q| < vr}. Note the change of the condition on the angle 
from Sin-i; the (double) critical point of (j){z; 1) is 2; = 1. As before, if z is in 
a compact subset of Sout,i, then there exists c > such that |e^'^(^'''')| < e"'^ 
uniformly in 7 G [1 — (5, 1 + 5]. 

Now define T by (80) as before. In doing so, we take Cm^i and Cin-i 
to lie in the annulus pi < \z\ < 1, and take Cout.i and Cout-i to lie in the 
annulus 1 < |z| < l/pi- Then the jump matrix in (81) satisfies 



(86) 



Vt{z) = I + O {e~^') 



uniformly for 7, 7 G [1 — (5, 1 + 5] and for z in all the contours except for 
{Cin-iUCout-i)<^{\arg{z)\ > vr-ai} and (Q„,i UCont,i) n{|arg(z)| < ai}. 



The parts of the contour where (86) is not valid are handled by introduc- 



ing local parametrix that can be solved by the RHP for the Painleve II 



equation (see Section 10). Such a 'Painleve parametrix' was introduced in 
the analysis of on a similar orthogonal polynomials but with a contin- 
uous weight. A drawback of the analysis of [4J was that the parametrix 
was solved asymptotically rather than exactly as in other cases such as 
|2H [22] . The exactly matching Painleve parametrix was constructed later 
in |19| . The construction of [19] requires, in the context of this paper, that 

a recent paper [11], a different approach 
to the exact construction of the Painleve parametrix was introduced. This 



construction has the advantage that it works for all 7 (and 7) in regime ( 72 ). 
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We seek a global parametrix in the form 
(87) A(z) = 



Ai(z) z&Ui 
A_i(z) zGU-i 
I elsewhere 



where U±i are sufficiently small, fixed size, neighborhoods of ±1. Later we 
will fix the size of U±i first and then choose 6 small enough so that U^i 
contains {Cin-i U Cout~i) H {|ar5((2;)| > tt — ai} and Ui contains {Cm^i U 



Cauti) n {|ar5(2;)| < ai} so that (86) is valid for all z in the contour of T 
except for in U±i. 

6.1. Local Models near 1 and —1. In order to construct exactly matching 
parametrices A-ti, we need to introduce Langer transformations which map 



the local phase functions and (j) to the Painleve phase (213) in lA-i and lAi 
respectively. 

The phase 6{z;^) is analytic in z in the neighborhood + 1| < 1 (and 
entire in 7) and admits the expansion, 

(88) e{z-^) = {l-j)(z + l) + ^{z + lf + ^^{z + lf + {{z + 1)') . 

At the critical value 7 = 1 the expansion degenerates to a cubic at leading 
order; for values of 7 near 1 the cubic unfolds either into three real or one 
real and two complex roots near z = —1. The double critical point-double 
root of 6'{z; l)-unfolds into a pair of simple critical points near z = —1: 

d f) 

(89) — = z± = --f± - 1 
dz 

Note that the relation (f){z; 7) = —9{—z; 7) implies that (p admits a similar 
expansion about z = 1 with the same structure. 



As the cubic coefficient in (88) is bounded away from zero (note that 
7<l + (5<3/2)we make use of a classical result of [16] to introduce new 
parameters 0(7) and 6(7) such that the relation 

(90) p(z-jf + a{j)f{z;j) + b{^) = -i9{z-n) zeU^i 

defines an invertible conformal mapping f = f(z) from a sufficiently small, 
7-independent, neighborhood U-i onto f{Ui) such that the parameters a 
and b depend continuously on 7 near 1. It was shown in [16] (see also [28]) 
that there exist 5i > and a 7-independent neighborhood ZY_i such that the 
above map is conformal in U-i for all 7 G [1 — (5i, 1 + 5i] if the critical points 
f± = lb -y/— a/2 of the left hand side, seen as a function of /, correspond to 



the critical points z± of 9{z;-y). This means that the left hand side of (90) 
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z^l 
f=fi/3f(z) \ 





Figure 7. Images of the contours near z = ±1 under 



evaluated at f = f± should equal to the right hand side of (90) evaluated 
at z = z±. These two conditions determine the parameters a and b as 



(91) 



b{-f) = ^[e{z+-n) + 9iz.;^)], 

(-«(7))'/' = |w^+;7)-^(^-;7)) 



Since 9{z+;'j) = —9{z-;j) = ^^7^ — 1 — 7 log(7 + -^7^ — 1), we have 

(92) 6(7) = 0. 

There are three choices of branch of 0(7). We choose the branch so that 

(93) a(^) = _^3i(^y^^_^ log (^^ + 7^31:^^]'/' 
satisfies the power series expansion 

(94) a(7) = 2(7 - 1) - ^(7 - 1)' + O ((7 - 1)') . 

To verify this, it is useful to note that ^ [\/7^ — 1 — 7 log(7 + \/7^ — 1)] 
— (7^ — With this choice of a, we have 

(95) /(.;„ = !<I^(. + l) + !(2_i) + 



+ 



1 /3-27 f7-l)' 



6i 



■) (z + l)3 + 0((^ + l)4). 



Inserting (94), we obtain 



(96) /(z;7) = -(z + l) 



l + liz + l) + ^iz + lf + 0{iz + lf) 



7 



+ ^(7-l)(^ + l) 



l + l{z + l) + 0{{z + lf 



+ 0((7-l)^(z + l)). 
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Define the rescaled coordinates (Langer coordinates) ( = ({z; 7) = t^^^f{z; 7) 
for z £ U-i, and set 



(97) 



.(7) = t'/V7). 



Then (see ( |2T3| )) 
(98) 



te{z;j) = i{-C + sC =iepn{C,s) 



zeU^ 



1- 



We note from (94) that for the parameters {n,m,t) in the regime (72), 
(99) 5(7) > -2L 

for all large enough t. We also have 

(2t2/3(^_l))^ 



(100) 5(7) = 2*2/^(7- 1) 



60 



t-2/3 + 0^2/3(^_l^ 



We introduce similar coordinates in Ui. This can be easily achieved by 
noting the symmetry (/'(z,7) = —6{—z,^). We set Ui = —U-i and define 
f{z;^) := —f{—z\^) for z G lAi. Then we find, with the same choice of o 
and 6, 

4 



(101) 



-/(z;7)-' + a(7)/(^;7) 



Z £ Ui. 



Defining Q = C,{z; 7) 
before, we obtain 

(102) t4>{z-^) = i 

Note the symmetry 
(103) 



= t^/^/(z;7), z G U-i, and s = 5(7) = t^/^a(7) as 
[C(^; 7)' + smiz; 7)1 = iOpiiiC, s), z^U.i. 



C{z-^) = -C(-^;7), 



z e Ui. 



We take 6 such that 5 < min{l/2, 5i} where 5i we introduced in defining 
/ in (90). Then consider the parameters {n,m,t) satisfying (72). 



Consider the image ofU-i under the map z 1— )■ C(z;7). From (96), we find 
that there exists 62 > such that for 7 E [1 — J21 1 + '^2], C(^-i; t) contains 
a disk centered at and of radius > O {t^^^) in the (^-plane. The same 
holds for Note that from ([96]), the image contours C(C'_i^j„/out) 

are oriented left-to-right and the image contours C(Ci,m/oMt) are oriented 
right-to-left as depicted in Figure [7] 

We now use (" to map the local contours and jump matrices inside U±i 
onto the jumps of the Painleve parametrix, RHP 10.1 We locally deform, if 
necessary, the contours C±i^in/out so that the image contours C(C±i,m/out 
li±i) become the rays Fj, i = 1,3,4,6 described in (206) and we extend 
C±i,m/ont 1^ (^-1 U^i) to the rest of C_|_i j„/o„j so that the estimate (86) 
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Figure 8. The jump contours for the residual R(z). The dashed 
lines represent contours on which the jumps are exponentially near 
identity. 



holds for z on the contour outside of U±i. The exact shape of the contours 
are not important. Reorienting the image contours, if necessary, to go from 
left-to-right and using (98) and (102) the image contours and jumps are, 
up to a conjugation by a constant matrix, exactly those of the Painleve 
parametrix, RHP 10.1[ 

Let ^'(Ci s) be the solution of the Painleve II model problem, RHP 10.1 
Set (72 = ( ^ "^^ ) and recall that 173 = ( q 1 ) • Taking into account the orien- 
tation of C,{C±i^in/out l^^iti)) we define the local models 

A_i(z) = A_i(z;7) :=a2a3"*(C(z;7);s(7))a3V2, z G 

Ai(z) = Ai(z;7) :=CT2*(C(2;7);s(7))f^2, z eUi. 

Note from the symmetries (214) and ( 103[ ) that these two models are related 



as 



(105) 



Ai (z, 7) = o-iCTgA-i (-2,7)0-3 cji, zeUi. 



From (98) and (102), A-i-i(z) satisfies the same jump condition as T(z) in 
U±i respectively. 

Define the ratio of the global parametrix to the exact problem T(z): 
(106) R(z) = T(z)A-i(z). 

Then R(z) has no jumps inside U±i but gains jumps on the positively ori- 
ented boundaries dU±i. Let S^j = U T^out U Cin,±i U Co«t,±i\(ZYi VJU-i). 
Then R satisfies the following problem 

Riemann-Hilbert Problem 6.1 for R(z): Find a 2 x 2 matrix R(z) 
such that 

1. R(z) is analytic in CXT^r where S^j = S^j U dUi U dU-i. 
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2. 'R{z) I as z ^ oo. 

3. The boundary values o/R satisfy the jump relation R+ = R_Vr where 

Ai(z)-i zedUi 

(107) Vr{z) = { A_i(z)-i z G dU-i 

Vt{z) zG^% 

The jumps of R(z) are now everywhere uniformly near identity. In fact, 
for the parameters {n,m,t) in the regime ( [72] ), it fohows from (86), 

(108) ||FR-/||i^(20^) = O(e-^*), 

and from ( |104 ) and (216) that (recall that C(^±ii t) contains a disk of radius 
>0{t^/^) for an7G [1-5,1 + 5]) 

(109) \\VR-I\\L^^u^^) = 0{t-'/^ 

(We will use a better estimate for the latter below.) The above estimates 
establishes that R falls into the class of small norm RHPs for any sufhciently 
large t. Let C_ : -L^(Sij) — L^(Sij) denote the usual Cauchy projection 
operator and define 

(110) CvMm:^c^lfM(VR-i)] = ^J /MlMii^d,, 

and 

(111) ^.im ^ / »M^d.., 

which maps / G L^{T,r) to an analytic function in C\Sjj. Then as C_ is a 
bounded operator whose operator norm is uniformly bounded (see e.g. 
^13j ) and the contours E/j are finite length it follows that ||C'yfl||j;^2_i.^2 = 
O (f~^/^) for large t which guarantees the existence of a unique solution 
to (1 — Cvji) fJ- = L Once the existence of fi{z) is established it follows 
immediately from the general theory of RHPs that 

(112) R(.):=/ + M;.](.) = /+-L/ l^HiVnH - D 

2vrz y„ w — z 



is the solution of RHP [6Tl 

Unfolding the series of transformations Yi— s-Qi— s-Si— T-T^-s-Rwe have 

r 



Y(0) = R(0) (^_^ and from ([67]) it follows that 

(113) vrn,m(0; t) = -Ri2(0; t, n, m) = R2i(0; t, n, m). 

We now evaluate R(0; t, n, m) explicitly for the first three terms in the as- 
ymptotic expansion. But we first consider the corresponding RHP for the 
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continuous weight in the next subsection. We will compare the discrete 
weight problem to the continuous weight problem. 

6.2. Analysis of the continuous weight problem. A streamlined ver- 
sion of the above procedure reducing the discrete problem, RHP |4.1[ to 
small-norm form can be used to study the continuous weight problem, RHP 
4.2 Using the same g-function used in the discrete case we define Y°° i— )• S°° 
)) replacing Q with Y°°. The new RHP for S°° features the single 
z;7) defined by (77) which we recall has a critical value at z = — 1. 



In the 'exponentially small regime' (71) estimate (78) holds and just as in 



Proposition 5.1 we have in the end 



(114) 7r„,oo(0;t) = O (e"^") for (n,t) satisfying ([71). 

In the Painleve regime ( |72[ ), by introducing a simplified version of the trans- 
formation (80), using only the factors appearing in Q±fl to open lenses, one 
defines a transformation S°° i-^ T°°. The problem for T°° is then approx- 
imated by a parametrix which is identity outside a neighborhood U-i of 
z = —1 and inside U-i is approximated by the same model as the discrete 
case, A_i(z) defined by (104). The result is a small norm problem R°° for 



_i(z) defined by (104). The result is a small norm problem R' 
the continuous case where 



(115) 
where, 

(116) 



1 
27ri 



fi°°{w)iVR^iw)-I) 



dw 



w 



Vroo(z) 



|A_i(z)-i zedU-i 

Moreover, the continuous weight orthogonal polynomial vr^(O) is given by 
(117) 



)(0;t) 



-Rr2(0;t,n) 



R^ (0; t, n) for (n, t) satisfying ( 72 ) 



6.3. Expansion of R(0). In this section we calculate the asymptotic ex- 
pansion of 

1 f ^i{w){VR{w) - 1) 



(118) 



R(0) =/ + /C/jM(0) =/ + 



-dw 



27ri ./^„ w 

up to order O {t~^). We begin by representing /x using its Nuemann series 
expansion. 



(119) 



k=l 



which, due to (108) and (109) , convergences uniformly and absolutely. In 



both (118) and (119) the dominant contribution to the integral comes from 



the boundaries dU±i. In fact, denoting by Pq the projection operator onto 
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^R\{dU-i^dUi), we find from ([TOS]) that ||CijPo|lL2(2^)^i2(s^) = O (e"^*) 

and ||/C^Po|lL2(2,)^^.(s,) = 0(e-'^*). 

Denoting by P±i the projection operator onto dhl± respectively, define 
C±i := CrP±i and /C±i := JCrP±i: for any / G ^^(Sk), 

J_ / /H(FflM- /) 

/■ /H(VrM_- ^) 



(120) 



C±i[f]iz) 



dw, ^ G Sij, 
dw, zf^YjR. 



Then we find 

(121) R(0) = /+(/C_i + /Ci)M(0) + O(e-^*) 



where 
(122) 



/x(^) = / + X;(C-i + Ci)'=[/](z) + 0(e 



k=\ 



Recall 5(7) defined in (97). Introduce the shorthands s = 5(7) and s 



5(7). Using ( |2T6| ), ( [217| ), and p04| ) we have 
(123a) 

y3(^) , ri ( e-'^ol'-l^'^^ 



Vr{z)-I= { 

with 

(123b) if2{s) 
¥'3(s) 

and, 
(123c) 



■^iCs) I 02(5) I 03(5) I 



(4/3 



z G dUi 



1 

2i 



1 



(2i)2 
1 

(2i)3 



in(s)2 - lg(,)2 (_l)n(5(,)^(,)_g'(s)) 



(-l)ng(s)«(.)-g'(.)) in(. 

a(s) (-l)"/3(s) 
-i-irm -a{s) 



,2 1 



0fc(s) = cT^Vfc(s)f^3", fc= 1,2,3, 



where g is defined by (210) and n, a, and /3 are defined in (216c)- (216e). 

It follows from inserting the above expansions into (121) and (122) that 
each iteration of Ci or C_i introduces a factor of t^^/^, thus we are led to 
an expansion of the form. 



N 



(124) 



R(o) = i + ^i?Wr^/=^ + o 



k=l 



t(A^+l)/3 
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where R^^^ := t^/'^{)Ci[I]{0) + /C_i[/](0)), 



(125) i?W:=t'=/3 J2 (/Ci +/C_i)C7^[I](0), k>2. 

re{-i,i}'=-i 

Here Cf is a multi-index understood as follows: given r = {ti,T2, . . . , r^) E 
{ — 1,1}'^ we define := Cr-i^Cr2 ■ • • Cr^.- Though we have suppressed the 
dependence, each R^^^ is a function of t. Moreover, since both s and the 
coefficients in the expansion (95) depend on 7, each R^^^ = O (1) with an 
expansion in powers of t~^^^. 

At each order we can split the composition of Cauchy integrals into three 
parts. Define 



(126) 



r: 



(fc) 

X 



Rik) _ rW _ R 



Note that from definition, R^^^ = 0. Intuitively, the first two "pure" terms 
contain the expansions of the continuous weight polynomials related to the 
marginal distributions while the last term contains the "cross" terms. This 
can be made concrete as follows. Let R±i(0) and Rx(0) denote the sum of 
each type of contributions to R(0): 



(127) 



Rp(0) := / + 



00 (fc) 

EJXp 

k=l 



P 



1,-1, X 



Clearly, Ri(0) and R_i(0) are the values at origin of normalized Riemann- 
Hilbert problems whose jump conditions are 



(128) 



(R_i)+(z) = (R_i)_(^)A_i(z,7)-i z G dU^i 
(Ri)+(z) = (Ri)_(z)Ai(z,7)-i zedUi 



Recalling (115) and (116) we see that R_i(2;) and R°°(z;t, n) have the same 



jump condition up to the exponentially small contributions from S/joo\(9ZY_|. 
Hence 



(129) 



R°°(0;t,n) = [l + 0{e-'')] R_i(0). 



Also from (105), the jump of Ri(2) is same as that of aia^1il°°{0,t,2m 
n)a'!^ai, and hence we find that 



(130) 



cJifT^R°°(0,t,2m-n)fTjai = [l + O (e"'^*)] Ri(0). 
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Therefore, from (117) it follows that 
(131) 



vr„,oo(0;t) = -(R_i)i2(0) + O {e-"') , 
7r2m-n,oo(0; t) = (-l)'^(Ri)i2(0) + O (e^^*) 



and hence from (113), (124) and (127), we find that 



(132) 7Tn,m{0; t) = 7r„,oo(0) - (-l)"vr2„_„,oo(0) - (Rx)i2(0) + O (e-^*) . 



(k) 

From ( 127), we now need to evaluate -Rp , p = —1, 1,X, k = 1, 2, 3. This 
calculation is a straightforward but lengthy application of residue calculus. 
We summarize the result of the calculations which follow directly from the 



definitions (126), (120), (123), making use of the expansions (95) and (100). 



It is helpful to note that the symmetry (105) between Ai and A_i implies 
that 



(133) 



where IC^J^'^^ and C^'^^'^ denote /C_i and C_i with 7 replaced by 7 respec- 
tively, and T is the operator defined by 



(134) 



Tfiz) ■.= a,a^f{-z)a^a,. 



In particular, note that TI = I, r[''^ = TR^^ll^^x^. 
Let Err and Err denote any terms satisfying 

(135) Err = O (^e~^o\s{y{T)r/^^ ^ ^ q (^^-co\s{^{r))\y^\^ _ 

Denoting by [A, B] and {A, 5} the commutator and anti-commutator of 
matrices A and B respectively, we find from an explicit evaluation that 



(making use of (217)) 



R 



(1) 



2i 1 



1 



(7-l))Ui(5) 



{2if 



(136a) 



R 



(1) 



3QV, -,,jr..-, 20t2/3 
(^|7 - Ip -h t"^''^|7 - 1| + t"^) Err, 

-2.ll-i(7-l)))*.W + ^*3W 



I7 - ir t-^i-'Yi - 1| + i~ Err 
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(136b) 



2 20*1/3 
+ (\l-l\+t-^^'^) Err, 



lOil/3 



20ti/3 

+ (I7- 1| Err, 



10*1/3 



(136c) 



^ M|)'^,(,)3 + f 1^ - 1| + t-V3) E,,^ 



R 



(3) 



20 

3(2i)3 

20 



(I7- 1| ^^^^^ 



(137a) 



(2) 

•X 



I7 - 1| + t~2/^) £^rr + (I7 - 1| + t"^/^) ^rr, 



(137b) = ^ {v=i(s)</>i(5)} {Ms) - Ms)) 



+ f|7 - 1| + M/^\ Err + f|7 - 1| + M/'^\ Err, 



Recall that i?^^ = 0. Note that 



(138) {Ms)Ai{s)} = 2{u{s)u{s) - {-irq{s)q{s))L 

From ( pTj ) and ( [l32| using ([123]) and p^-jpT]), we obtain 
Proposition 6.1. Set 



(139) 



9i{y,y) ■= -j^{u {y)q{y) + u{y)q {y)) , 
g2(.y,y) ■= ^iq{y)u' (y) + q' {y)u{y)) . 



Let T^n,m{z) he the orthogonal polynomial given in (67). Let T^n,oo{z) he the 



orthogonal polynomial given in (69). There exists 5 > ^ such that for any 
fixed L > 0, if 



(140) 2t-Lti/3 <n<2t(l + (5), 2t - Lt^/^ < 2m - n < 2t{l + 6), 
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then there exists constants cq > and to > such that 
(141) 

7rn,m(0;t) =TTn,oo{0;t) - (-l)"7r2m-n,oo (0; t) 

^ g,{s{j),sm-{-irg2{s{j),sm 

/or a// 1 > to, where 

(142) " , 2m-n 



7 := 



n 

2t' 



2t 



and is defined in (97) which satisfies (see (lOO)J 

(143) S(^) = 2t'/'(^ - 1) - i^^!^!^ilLlll!t-2/3 + (^i2/3(^ _ -^)3^ _ 

We also have 
Proposition 6.2. For t > to, 



(144) 

where 
(145) 



(-l)"7r„,o,(0; t) =^q{s{j)) (l - ^] + h{s{^)) 



+ O ((t-4/3 + i-2/3|^ _ l|)g-c„K7)r j 



h{y) := ^u{y)q'{y) - ^q^ - ^Vqiv)- 



7. Proof of Theorem 11.11 and Corollary 11.11 

We now evaluate the asymptotics of P{CRf < k, NE^ < j} when 

(146) j = [t + 2-^xt^/^], k=[t + 2~'^x't''^/% 

where x,x' G M are fixed, and [a] denotes the largest integer no larger than 
a. We define xt and x[ by 

fl47) X. - + - x' - + - 

SO that 

(148) 2i + 1 = 2t + 2k + I = 2t + x[t^/^ . 

Then xt = x + 0{t-^/'^) and = x' + ©(r^/^)^ 
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From Proposition we have 



(149) logP{CRt < A;,NEt < j} 



/ vr2,+i,^(0; T)dT+ / / Qf{T)dTds 

Jo Jo Jo 



where 



(150) 



and 



Q-(r) = -7^-(r) -5f (r) +7^-(r)5f(r), 



(151) 7^™(r) := vrs,- ^(0; T)7r2,-+2,m(0; r), 5f (r) := |7r2j+i,„,(0; t)|2. 



From Proposition 5.1 (substituting r for t in (84)), we find that the above 
integrals away from the interval [(1 — €)t,t], for any fixed e > 0, are expo- 
nentially small in t: 



(152) log P {CRt < fc, NEt < j} 

= / 7r2,+i,m(0;r)dT+ / / Q^{r) dr ds + 0{t 

Jt(l-e) Jt(l-t) Jt(l-t) 



-ct\ 



We can take e > small enough so that Proposition 6.1 is applicable to 
7^2j+i,m{^;T) for £ = 0, 1,2 and r G [(1 - e)t,t]. 
Now by the same argument, we have 



(153) logP{NEt<i} 

/i^ /i^ fig 

= / vr2j+i,oo(0; r)dr+ / / Qf (t) dr ds + 0{e 

Jt(l-e) Jt(l-e) Jt(l-e) 



and 



(154) logP{CRi<A;} 

vr2fc+i,oo(0; T)dT+ / [ Q^[r)dTds + 0{t 

Jt(l-e) Jt(l-e) 



f 

Jt{l-t) 



Consider 



(155) logP{CRt < A:,NEt < j} - logP{NEj < j} - logP{CR < A;} . 
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We first consider the three single integrals. From ( 141 ) applied to n 
2j + 1 and t replaced by r, we have 



(156) 



/ [vr2j+i.m(0;T) - 7r2j+i,oo(0;r) - 7r2fc+i,oo(0; r)] dr 

Jt(l-e) 

(^-4/3 + r-2/3|^(^) _ l|)e--o|.(7W)l^/^dr^ 



1 



t(i-e) r 
+ 



where 
(157) 



7(r) : 



2i + l 



7(t) : = 



2k + I 



2r ' " ' 2r 
Changing the integration variable r i— )• r/ as 

(158) T = t- 2-^r]t^/^. 



the integral involving gi in (156) becomes 



(159) 



1 



2*2/3 



5ri(s(7(r),s(7(r))) 



dr] 



l-2-ir/t-2/3- 



Note that from ([100|, 
(160) 



s(7(r)) = (x; + r?) + 0(r?V2/3). 

Also note that from its definition, gi (xq+t], Xq+t]) is integrable for r] G [0, oo) 
for any fixed xo,Xq £ M. Thus, we obtain that the integral (159) equals 

1 



(161) 



2t2/3 J, 



9i(xt 



+ r],x[ + 7])drj + (t-^/^) . 



The integral involving g2 in (156) equals the same integral with gi replaced 
by 52- On the other hand, it is easy to see that the error term in ( |156 ) is 

(162) O (t'/^ r t-^'/^il + \xt + T]\)e-'^^^'+^^'^'dv) = O (t-i) . 



Thus, replacing xt and x[ by x and x' , which incurs an error of order 
O ( fl56l ) equals 



(163) 



1 f°° 

2^2/3 [gi{x + 'n,x' + ri)+g2{x + r],x' + 'n)]dr] + 0{t'^) . 



Now inserting the definition (139), we can perform the integration and we 
find tl 

(164) 



find that (156) equals 
-1 



4i2/3 



[u{x)q{x') + q{x)u{x')] + O {t~^) . 



CROSSING AND NESTING 



39 



We now consider the part of (155) that come from the three double inte- 
grals. We need to evaluate Q^(t) - Qf(r) - Q^(r). Setting 



(165) 



± 2j + 1 ± 1 1 

' ^ ^ 2r ^ 2r' 



we see from (|100|) that 
(166) 
Let us set 

(167) ^:=s(7(r)), i := s{^{t)) 



.(7±(r)) = .(7(r)) ± -^^ + 0(t-i/^(7(r) - 1)). 



to ease the notational burden. Then, (144) implies, using (217), that 

(168) ^2i+i±i,oo(0;r) = - 7r2j+i,oo(0; r) ± q'i^:^ + \^"(0^ 

+ T~^^^ Error. 

where throughout the rest of this section we use the notation Error to 
denote any term satisfying 

Error =0 ((l + t^'^'Ht) - l|)e-^ol^(^M)l'^'' 

(169) ^ 

+ 0((l + r^/3|7(r)-l|)e-'=»l^(^M)l''' 

Note that 
(170) 

't{l-e)Jt(l-e) 

Also, note that from (144), (168) implies, in particular, that 

(171) vr2j+i±i,oo(0; r) =q{^) + t'^'^ Error, 

and clearly the asymptotics (168) and ( 171| ) also hold when j is replaced by 
k and ^ is replaced by ^. 



I [ Error drds. = 0{t^/-'), 

Jt(l-t) Jt(l-e) 



From (141), 



(172) 



k2j+i,m(0; t)P - |7r2j+i,oo(0; r)p - |vr2A;+i,oo(0; r)|^ 

= 27r2j+i,oo (0; r )7r2fc+i,oo (0; r) 
2 

+ - 

T 



[51 (C, + 92{i, 0] [7r2i+i,oo(0; r) + 7r2fc+i,oo(0; r)] 
+ T~^/'^ Error. 
Thus, from ([l44]), 

5™(r) - 5f (r) - 5^ (r) = 27r2,-+i,oo(0; r)^2fc+i,oo(0; r) 
(173) 2 



r4/3 



[91 (e, + 52(e, I)] [g(e) + g(l)] + r-^/^i?' 



rror. 
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Similarly, using put and p7T| ), we obtain 
7^f(r)-7^f(r)-7^r(T) 

= -'^2j,oo (0; r )7r2fc,oo (0; r) - 7r2j+2,oo (0; r )7r2fc+2,oo (0; r) 



(174) 

and 
(175) 



+ :;473 [aiii. 1) - 52(e, 0] W) - q{i)] + t-'/'e^ 



n^{r)sf{T) - nf{T)sf^{r) - nr{T)s^{r) 

-2 



rror. 



- ^4/3 9(09(6 +r ^''^ Error. 

Therefore, since 

7r2j,oo (0; r )7r2A;,oo (0; r) + 7r2j+2,oo (0; r )T^2k+2,oo (0; r ) 
(^;^7g-) - 2vr2j+i,oo(0; r)7r2fc+i,oo(0; r) 

= [9(e)g"(6 + q"{0<l{b + 2g'(09'(0] + r-^'I^Error, 

we obtain, by using the definition of gi , g2 and by using the fact that 
and 2(7(7' = that 

(177) Q-(t) - Qf (r) - Q^(r) = ^W(e, f) + r-^/^i^rror 



where ^ := s(7(r)), ^ := s{^{t)) are defined in (167) and we have set 
m, i) ■=n"iCMO + 2n'(0^'(l) + uiOn'ii) 



We insert (177) into the integral 

(179) r r [QT{r)-QT{T)-Qr{r)]dTds, 

Jt{l-t) Jt(l-e) 

and evaluate it by changing variables t rj, t = t — l^^r/t^/^ and s i— )• 
s = t — 2~^(^t"^/^, as was done for the single ingtegrals. Noting that 

(180) U{C + r],C + r]) := ^ [ui^ + r/)u(e + v) + 9(C + ^)'?(C + r])] , 



the integral can be evaluated and we find that (|179|) equals 
(181) 



1 



^^2/3 Hx)u{x') + (z(x)g(x')] + O (t-1) . 

The error term O {t~^) follows from ( 170| ). 

Combining (164) and (181), we obtain 
(182) 



log 



P{CRt < X, NEt < x'} 
IP{CRt < x}P{NEi < x'} 



[q{x) — u{x)] [q{x') — u{x')] 

4t273 



+ o(t-'). 
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This completes the proof of Theorem We note that here the error term 
is uniform for in a compact subset of M (actually in any semi- infinite 
interval [xq, oo).) 

t-y^ + o(t-^). 



Corollary 



1.1 



follows if we show that Cov(CRt,NEt) 
This is obtained from Theorem |1.1| by using the dominated convergence 
theorem if we have tail estimates of P{CRt < x, NEt < x'} - P{CRt < 
x}P{N'Et < x'} as — )• oo since J^^xdF'{x) = —1. The tail as 
X, x' —7- +00 can be obtained from the analysis of this paper. For the other 
limits, we need an extension of the analysis of this paper but we skip the 
details in this paper. See [Ull] for a similar question about the convergence 
of moments using Toeplitz determinant. 

8. Proof of Theorems 11.21 and 11.31 
Here we evaluate the asymptotics of the marginal distributions P{CRt < 



j} for j as given by (146). We reuse as much as possible the calculations 
in the previous section. Note that by symmetry we have P{NEj < j} = 
P{CRt < j}. In the process of computing the marginal we will compute as 
a by-product asymptotics for F{Lt < i} along the way. 
Our starting point is to introduce the change of variables 

(183) T = t-2-^{r]-xt)t^/''^ s = t-2-^{C-xt)t^/^ 



into (153) where, as in the previous section, xt is given by (148). Note that 



this change of variables differs from ( 158| ) by a shift. Making the substitution 
we have, with j and k defined by (146) (recall (148)), 

log P {NEt < j} = Xi + X2 + O (e-^*) , 



(184) 
where 

(185) 



logP{Lt < 2j + 1} = 2X1 + O (e-^*) , 



Xt -^C 

Xt+2et^/^ 

1T2j + l,oo{0;T) dr]. 

Xt 



Qf{T)dr^dC, 



From ( 63 ) , there is an analogous formula for log P {Lt < 2j}, and the analysis 



below applies to this case too without much changes. We skip the details 
for this case. 

In order to compute expansions of the above integrals we need more de- 



tailed calculations than the previous section. Inserting (183) into (157), we 
have 



(186) 



7(r) 



i + lvt 



-2/3 



■qxt 
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Then (143), with t replaced by r, becomes 



(187) 



s(7(t)) =7]+ {—r] 



3 2 1_ W-2/3 + ( 3^-4/3 



Inserting these into ( 144 ) we have 



(188) -7r2j+i,oo(0;r) = 



Hv) + ( ~ ) (^('^^ + w'(^)) 



+ ( r'^/^ Error 



and it fohows from (150), (151) (when m = oo), and (a shght improvement 



of) (168) that 



(189) Q-(r) = -2t-2/3g(r?)2 



-4/3 



4q{rj)h{r]) + ( - '^Xt ) {m {r])q{Ti) + ^(77)^) 



j{r])q"{r]) - q{r]) - q{r]) 



Here h is as given in (145). In both the above formulae the Error term is 



as defined in ( 169 ) and we recall that its integral introduces terms of order 



O (^2/3). Now using the identity = n + {q') — rjq'^ and using the fact that 
q"^ = u' , 2qq' = u" , and q" = rjq + 2q^ it is direct to check that the terms in 



square brackets in (188) and (189) can be expressed as perfect derivatives. 
We find that 



(190) 



-^2j+i,oo(0;r) = ^q{r]) + hlM + O {t-^/^Error^ , 



2 



u{7]) - 4nU2i7]) + O ir^/'^ Error 



where 



(191) 



5 di] 
1 d^ 

5 dr]'^ 



u{r])q{r]) - q'(r]) + — {9rj - Wxt) r]q{r]) 



12 



uiv)"^ - q{f]f + 77 (9r/ - lOxt) ??u(??) 
6 
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Inserting this formula into ( 185 ) and ( 184 ), we obtain with x^*-* and xt defined 
by (26) and (13) respectively, 



(192) logP {it <2t + t^^^x} = XogFcuEi^ 



10*2/3 



x(x)^ — q{x)'^ x^u(x) 

6 



+ o(t 



-1^ 



and 

(193) log P {nE^ < t + 2-h^/^x} = log F{xt) + ^ + O {r^) 
where E = E{x) equals 



(194) E :-- 



20 



(u(x) - qix))"^ + 2 iu'ix) - qix)) + iuix) - qix)) 

o 



It is easy to check that 2{)E{x)F{x) = -4F"(x) - WF'{x) and {u{xf - 
q{xY -\x'^u{x))Fgue{x) = F2f;^(x) + |x2F^^^(x)g Theorems 
follow immediately. 



1.2 



and 



1.3 



9. Proof of Corollary 11.21 
For a sequence {a„}^0) consider its Poissonization 



(195) 



2\n 



it') 



-an. 



n=0 



A de-Poissonization lemma is that if (a) < < 1 and (b) a„+i < a„ for 
all n, then we have for s > 1 and n > 2, 



(196) 

where, 

(197) 



1 



< an < 



1 



Hn--=n + 2i/snlogn, i/„ 



n 



2\J sn log 



n. 



Lemma 2.5 of |33j is stated for the case when s = 2, but the proof can be 
modified in a straightforward way to obtain the above estimates. 

The de-Poissonization lemma can be applied to := P {cr„ < /c, ne„ < j} 
due to the following lemma. 

Lemma 9.1. For each n > 0, and k,j > 0, 

(198) P{cr„+i < k, ne„+i < j} < P{cr„ < k, ne„ < j}. 



We would like to thank Craig Tracy for pointing out these relations. Relations like 
these and many others can be found in ,47, . 
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Proof. Since P{cr„ < k, ne„ < j} 
(199) gk,j{n) := #{Af G A1„ : cr„(M) < A;,ne„(M) < j} , 



we need to show that gk,j{n + 1) < (2n + l)5fcj(n). The set A^n+i of 
complete matchings of [2{n + 1)] is the union of (2n + 1) disjoint subsets 
1, . . . , 2n + 1, where -M^ , is the set of complete matchings 



M 



n+l' 



of [2(n + 1)] such that 1 is paired with £ (i.e.(l,£) is an element of the 
matching). By removing the two vertices 1 and £, and then relabelling the 
vertices, there is a trivial bijection : A^^+i ^ M.n- Clearly, cr„+i(M) > 
cr„(/£(M)) and ne„+i(M) > ne„(/^(M)) for M £ M^+i- This implies that 
gk,j{n + 1) < (2n + l)gkj{n). □ 



Hence, since (see (16)) 



'{CRt<k,NEt<j} = e~''/^Yl 



(200) 



n=0 



n=0 



'{cr„ < k,nen < j} 
'{cr„ < k} , 



we find that for each s > 1, n > 2, and j, k > 0, 
(201) 

F {cvn < k, ne„ < i} - P {cr„ <k}F {ne„ < j} 

< IP {cRv2i^ < k, NE^ < j} - P {CR^ < A:} P {nE^ < j} + 4n-^ 

When k = V2n + 2-ix(2n)i/6 and j = V2n + 2-^x'{2ny/^, from Theo- 
rem [lT| the right-hand-side of (201) is less than or equal to 

(202) 



CR 



< k 



}r{ 



NE 



< j 



+ 4n-' + o( n-^/^ 



Now we use Theorem 1.2 to estimate each of the above probabilities. Note 
that 



(203) 



2n-h2-ix(2n)i/6 - ^2^ 



4 v/sn log n / v/log n 

2; H ^ TTTT 1- C 



2-i(2z^„)i/6 (2n)i/6 ni/2 

When i^n is replaced by then the first plus sign on the right-hand-side 
is changed to the minus sign. From this, it follows that (202) is bounded 
above by o(^) 



+ An ' . The lower bound is similar. Thus we obtain 



Corollary 1.2 
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10. A MODEL RHP: Painleve II 
Consider the coupled pair of differential equations for 2x2 matrix ^'(C, s): 



(204a) 
(204b) 



d^- 
d^- 

^d7 



(4C2 + s)[a3,^] + 



2q' 
AiCq + 2r 



4:iC,q — 2r 
-2q^ 



where as denotes the Pauli matrix and [*, *] is the commutator: 

[A, B] = AB — BA. The compatibility condition for this overdetermined 
system is that q = q{s) satisfy Painleve II q" = sq + 2q^ and r = q'{s). This 
is a representation of the Lax-pair for Painleve II equation introduced by 
Flaschka and Newell 1261. 




Figure 9. The contours Tj and regions Sj defining 5'((',s) 



Any solution of (204a) is an entire function of C- Let Sj,j 
denote the sectors 



1,...,6 



(205) 



Sj = ->C€C : < arg(C) < ^^vr 



6 



6 



and let F,- denote the outwardly oriented boundary rays: 



(206) 



2i-l 



- >,CeC : arg(C) - — 



-vr 



There exist a unique solution of ( |204a ) such that 
(207) = / + O (C"^) asC^ooinSj, 
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and constants a, , j = 1 , . . . , 6 such that for ( E F,- 



(208) 



^i+i(C) = ^i(C)L ' ^ 1, J even. 



Additionally, the constants aj satisfy 

(209) aj+3 = aj, 010203 + 01 + 02 + 03 = 0. 

The parameters Oj depend parametrically on s, q and r; in [26] Flaschka and 
Newell showed that the isomonodromic deformations, i.e., the variations of 
these parameters that keep the Stokes multipliers aj constant, are given by 
solutions of the Painleve II equation: q"{s) = sq + 2q{s)^ and r{s) = q'{s). 
Our particular interest is in the Hastings-McLeod solution of Painleve II 
which is the unique solution such that 



(210) g(s) = Ai(s) (1 + o (1)) as s — )• 00, q{^)'^\I~^ as s — )• — 00 




Let \I'(C; s) be the solution of (204a) with parameters s, q = q{s) and 
r = q'{s), where q{s) is the Hastings-McLeod solution, and let V denote 
the set of poles of q (of which there are infinitely many). Then ^'((^, s) is 
defined and analytic for C £ C\(Ci U C2) and s G C\P. It is known that 
there are no poles of q on the real line |32j . The Stokes multiplier for the 
Hastings-McLeod solution are 

(211) oi = 1, 02 = 0, 03 = -1. 

If we reverse the orientation of r3 and r4 and define Ci = Li U and 
C2 = r4 U Le, then ^'(C; s) solves the following RHP. 




Figure 10. The contours defining RHP 10.1 related to the 
Hastings-McLeod solution of Painleve //. The contours can be de- 
formed to the dashed lines without changing the problem statement. 
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Riemann-Hilbert Problem 10.1 (PII model RHP) Find a 2x2 matrix 
^'(C; s) with the following properties. 

1. ^'(C; s) is an analytic function of ( for ( G C\(Ci U C2). 

2. ^{C]s) = 1 + (C"^) as C ^ 00 and bounded as C ^0. 

3. The boundary values ^±{C,]s) satisfy the jump conditions 



(212) 

where, 
(213) 



*+(C;s) = vI/„(C;5) 
*+(C;s) = ^_(C;s) 



1 

1 



CGC2. 



epii = epn{c,s) = -c^ + sc 



We make two observations which we will need later. First, the symmetries 
— Ci = C2 and 9pij{—(, s) = —9pjf{(,s) imply that the solution ^'(C,s) of 
RHP |10.1 satisfies the symmetry 

(214) ^{-C,s) = ai^{C,s)ai, :=(?!). 

The second fact is that ^ admits a uniformly expansion in the limit as 
C — )• 00 as described in [23] . Specifically, we have 



(215) 



^{C;s) = I + 



C 



+ 



The error term O (C~^) here depends on s. For our purpose, we need the 
dependence on s for s bounded below. An analysis similar to Section 6 of 
|23j shows that given sq > 0, there exists a constant cq > such that 



(216a) 



*(C;s) 



-co|s| 



3/2 ' 



c 



The moments ipj{s) can be calculated recursively from inserting the expan- 
sion into (204b). The first three moments are 

Ms) 

(216b) Ms) 

Ms) 



1 

2i 



1 

W 
1 



u{s) q{s) 
q{s) u{s) 



2n(s) 2 
q{s)u{s) - q'{s 



a{s) 
Pis) 



Hs)' 

1' 

-Pis)] ^ 
-a{s) ' 



/(s)u(s) - q'{s) 



hi^)' 
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where 



(216c) u{s)= rqiO'd^, 

J oo 



(216d) = - "Jf + logF(,)' - IjKf^K. 



/3(s) = q'{s)u{s) - q{s)(^s + 



g(s)^ n(s)^ 



(216e) i3(s)=q'{s)uis)-qis){s+^+ ^ 

We note that the asymptotic analysis of the RHP for the Painleve equation 
imphes that for a given sq > 0, 

(217) ^,(s) = 0(e-^«W'''), i = 1,2,3, 



where cq can be taken as the same constant in the error term of (216a 



References 

[1] Adler, M., Ferrari, P. L., and van Moerbeke, P. Non-intersecting random 
walks in the neighborhood of a symmetric tacnode. ArXiv e-prints (July 2010), 
1007.1163. 

[2] Baik, J., Borodin, A., Deift, P., and Suidan, T. A model for the bus system in 
Cuernavaca (Mexico). J. Phys. A 39, 28 (2006), 8965-8975. 

[3] Baik, J., Buckingham, R., and DiFranco, J. Asymptotics of Tracy- Widom dis- 
tributions and the total integral of a Painleve II function. Comm. Math. Phys. 280, 
2 (2008), 463-497. 

[4] Baik, J., Deift, P., and Johansson, K. On the distribution of the length of the 
longest increasing subsequence of random permutations. J. Amer. Math. Soc. 12, 4 
(1999), 1119-1178. 

[5] Baik, J., Deift, P., and Rains, E. A Fredholm determinant identity and the 
convergence of moments for random Young tableaux. Comm. Math. Phys. 223, 3 
(2001), 627-672. 

[6] Baik, J., Kriecherbauer, T., McLaughlin, K. T.-R., and Miller, P. D. Dis- 
crete orthogonal polynomials, vol. 164 of Annals of Mathematics Studies. Princeton 
University Press, Princeton, NJ, 2007. Asymptotics and applications. 
[7] Baik, J., and Liu, Z. On the width of non-intersecting processes, in preparation. 
[8] Baik, J., and Rains, E. M. Algebraic aspects of increasing subsequences. Duke 

Math. J. 109, 1 (2001), 1-65. 
[9] Baik, J., and Rains, E. M. The asymptotics of monotone subsequences of involu- 
tions. Duke Math. J. 109, 2 (2001), 205-281. 
[10] Basor, E. L., and Ehrhardt, T. Determinant computations for some classes of 

Toeplitz-Hankel matrices. Oper. Matrices 3, 2 (2009), 167-186. 
[11] BORNEMANN, F. Asymptotic independence of the extreme eigenvalues of Gaussian 

unitary ensemble. J. Math. Phys. 51, 2 (2010), 023514, 8. 
[12] BORNEMANN, F. On the numerical evaluation of distributions in random matrix 

theory: a review. Markov Process. Related Fields 16, 4 (2010), 803-866. 
[13] Bottcher, a., and Karlovich, Y. I. Carleson curves, Muckenhoupt weights, and 
Toeplitz operators, vol. 154 of Progress in Mathematics. Birkhauser Verlag, Basel, 
1997. 



CROSSING AND NESTING 



49 



[14] Buckingham, R. J., and Miller, P. D. The sine-Gordon equation in the semiclas- 
sical limit: critical behavior near a separatrix. ArXiv e-prints (June 2011), 1106.5716. 

[15] Chen, W. Y. C, Deng, E. Y. P., Du, R. R. X., Stanley, R. P., and Yan, C. H. 
Crossings and nestings of matchings and partitions. Trans. Amer. Math. Soc. 359, 4 
(2007), 1555-1575 (electronic). 

[16] Chester, C, Friedman, B., and Ursell, F. An extension of the method of steep- 
est descents. Proc. Cambridge Philos. Soc. 53 (1957), 599-611. 

[17] Choup, L. N. Edgeworth expansion of the largest eigenvalue distribution function of 
CUE and LUE. Int. Math. Res. Not. (2006), Art. ID 61049, 32. 

[18] Choup, L. N. Edgeworth expansion of the largest eigenvalue distribution function of 
Gaussian unitary ensemble revisited. J. Math. Phys. 49, 3 (2008), 033508, 16. 

[19] Claeys, T., and Kulilaars, A. B. J. Universality of the double scahng limit in 
random matrix models. Comm. Pure Appl. Math. 59, 11 (2006), 1573-1603. 

[20] CORWIN, I., Quastel, J., and Remenik, D. Continuum statistics of the Airy2 
process. ArXiv e-prints (June 2011), 1106.2717. 

[21] Deift, p., Kriecherbauer, T., McLaughlin, K. T.-R., Venakides, S., and 
Zhou, X. Strong asymptotics of orthogonal polynomials with respect to exponential 
weights. Comm. Pure Appl. Math. 52, 12 (1999), 1491 1552. 

[22] Deift, P., Kriecherbauer, T., McLaughlin, K. T.-R., Venakides, S., and 
Zhou, X. Uniform asymptotics for polynomials orthogonal with respect to varying 
exponential weights and applications to universality questions in random matrix the- 
ory Comm. Pure Appl. Math. 52, 11 (1999), 1335-1425. 

[23] Deift, P. A., and Zhou, X. Asymptotics for the Painleve II equation. Comm. Pure 
Appl. Math. 48, 3 (1995), 277-337. 

[24] El Karoui, N. A rate of convergence result for the largest eigenvalue of complex 
white Wishart matrices. Ann. Probab. 34, 6 (2006), 2077-2117. 

[25] Ferrari, P. L., and Frings, R. Finite time corrections in KPZ growth models. J. 
Stat. Phys. 144, 6 (2011), 1123-1150. 

[26] Flaschka, H., and Newell, A. C. Monodromy- and spectrum-preserving defor- 
mations. I. Comm. Math. Phys. 76, 1 (1980), 65-116. 

[27] FOKAS, A. S., It-s, a. R., and Kitaev, A. V. The isomonodromy approach to 
matrix models in 2D quantum gravity. Comm. Math. Phys. 147, 2 (1992), 395-430. 

[28] Friedman, B. Stationary phase with neighboring critical points. J. Soc. Indust. Appl. 
Math. 7 (1959), 280-289. 

[29] Gessel, I. M. Symmetric functions and P-recursiveness. J. Combin. Theory Ser. A 
53, 2 (1990), 257-285. 

[30] GOLiNSKii, L. B. Schur flows and orthogonal polynomials on the unit circle. Mat. Sb. 
197, 8 (2006), 41-62. 

[31] Grabiner., D. J. Brownian motion in a Wcyl chamber, non-colliding particles, and 

random matrices. Ann. Inst. H. Poincare Probab. Statist. 35, 2 (1999), 177-204. 
[32] Hastings, S. P., and McLeod, J. B. A boundary value problem associated with 

the second Painleve transcendent and the Korteweg-de Vries equation. Arch. Rational 

Mech. Anal. 73, 1 (1980), 31-51. 
[33] Johansson, K. The longest increasing subsequence in a random permutation and a 

unitary random matrix model. Math. Res. Lett. 5, 1-2 (1998), 63-82. 
[34] Johansson, K. Non-intersecting paths, random tilings and random matrices. Probab. 

Theory Related Fields 123, 2 (2002), 225-280. 



50 



JINHO BAIK AND ROBERT JENKINS 



[35] Johansson, K. Discrete polynuclear growth and determinantal processes. Comm. 

Math. Phys. 242, 1-2 (2003), 277-329. 
[36] Johansson, K. The arctic circle boundary and the Airy process. Ann. Probab. 33, 1 

(2005), 1-30. 

[37] Johnstone, I. M. Multivariate analysis and Jacobi ensembles: largest eigenvalue, 

Tracy-Widom limits and rates of convergence. Ann. Statist. 36, 6 (2008), 2638 2716. 

[38] Johnstone, I. M., and Ma, Z. Fast approach to the Tracy-Widom law at the edge 
of GOE and GUE. ArXiv e-prints (Oct. 2011), 1110.0108. 

[39] Kamvissis, S., McLaughlin, K. D. T.-R., and Miller, P. D. Semidassical soli- 
ton ensembles for the focusing nonlinear Schrddinger equation, vol. 154 of Annals of 
Mathematics Studies. Princeton University Press, Princeton, NJ, 2003. 

[40] Karlin, S., and McGregor, J. Coincidence probabilities. Pacific J. Math. 9 
(1959), 1141-1164. 

[41] KONIG, W., O'Connell, N., and Roch, S. Non-colliding random walks, tandem 
queues, and discrete orthogonal polynomial ensembles. Electron. J. Probab. 7 (2002), 
no. 5, 24 pp. (electronic). 

[42] Krattenthaler, C. Growth diagrams, and increasing and decreasing chains in fill- 
ings of Ferrers shapes. Adv. in Appl. Math. 37, 3 (2006), 404-431. 

[43] Ma, Z. Accuracy of the Tracy-Widom limits for the extreme eigenvalues in white 
Wishart matrices. ArXiv e-prints (Mar. 2012), 1203.0839. 

[44] Moreno Flores, G., Quastel, J., and Remenik, D. Endpoint distribution of 
directed polymers in 1-|-1 dimensions. ArXiv e-prints (June 2011), 1106.2716. 

[45] Nenciu, I. Lax pairs for the Ablowitz-Ladik system via orthogonal polynomials on 
the unit circle. Int. Math. Res. Not, 11 (2005), 647-686. 

[46] Rains, E. M. Increasing subsequences and the classical groups. Electron. J. Combin. 
5 (1998), Research Paper 12, 9 pp. (electronic). 

[47] Shinault, G., and Tracy, C. A. Asymptotics for the covariance of the Airy2 
process. J. Stat. Phys. 143, 1 (2011), 60-71. 

[48] Stanley, R. P. Increasing and decreasing subsequences and their variants. In Inter- 
national Congress of Mathematicians. Vol. I. Eur. Math. Soc, Ziirich, 2007, pp. 545- 
579. 

[49] SzEGO, G. Orthogonal polynomials, fourth ed. American Mathematical Society, Prov- 
idence, R.I., 1975. American Mathematical Society, Colloquium Publications, Vol. 
XXIII. 

[50] Tracy, C. A., and Widom, H. On orthogonal and symplectic matrix ensembles. 
Comm. Math. Phys. 177, 3 (1996), 727-754. 

Department of Mathematics, University of Michigan, Ann Arbor, MI, 48109, 
USA 

E-mail address: baik@uinich.edu 

Department of Mathematics, University of Michigan, Ann Arbor, MI, 48109, 
USA 

E-mail address: rmjenkin@umich.edu 



